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Epigraph 



REALITY, n. The dream of a mad philosopher. That 
which would remain in the cupel if one should assay a phan- 
tom. The nucleus of a vacuum. 

-from "The Devil's Dictionary" , by Ambrose Bierce (1911) 
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Abstract 



We discuss recent progress in the understanding of the vacuum structure (effective super- 
potentials) of confining gauge theories with Af = 1 supersymmetry, in particular the- 
ories with softly broken Af = 2 supersymmetry. We show how the new techniques 
improve upon older calculations in non-supersymmetric quantum field theories. A com- 
mon feature of both approaches is that appropriate perturbative field theory calculations 
(e.g. using the background field method) give non-perturbative information about the 
vacuum structure of the theory. However, in supersymmetric theories, these results are 
often exact. 

The geometric engineering of supersymmetric gauge theories in string theory pro- 
vides powerful tools for studying gauge theories. Central to the analysis is a particular 
class of hyperelliptic curve, which emerges from the Calabi- Yau geometry of the string 
theory background and encodes the gauge theory effective superpotential. These curves 
may be rederived using other techniques based on zero-dimensional matrix integrals, the 
dynamics of integrable systems and the factorization of Seiberg-Witten curves, and we 
describe in detail how each technique highlights complementary aspects of the gauge 
theory. 

We find that the use of the spectral curve requires the introduction of additional fun- 
damental matter fields, which act as regulators for the UV divergences of the calculation 
by embedding the gauge theory in a UV-finite theory. Theories with < Nf < 2N C 

viii 



fundamental multiplets may thus be treated uniformly. We focus in detail on maximally- 
confining vacua of Af = 1 gauge theories with fundamental matter, and of gauge theories 
with SO and Sp gauge groups. Both cases require refinements to the basic techniques 
used for SU gauge theory without fundamental matter. 

We derive explicit general formulae for the effective superpotentials of J\f = 1 the- 
ories with fundamental matter and arbitrary tree-level superpotential, which reproduce 
known results in special cases. The problem of factorizing the Seiberg-Witten curve for 
Af = 2 gauge theories with fundamental matter is also solved and used to rederive the 
corresponding J\f = 1 effective superpotential. 
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Chapter 1 



Introduction 

A central problem in theoretical particle physics is to understand the nature of the strong 
nuclear interactions at low energies. A quantitative theory of the strong nuclear force 
(quantum chromodynamics, or QCD) has been known for over 30 years, but com- 
putational difficulties prevent accurate analytical calculations at low energies or long 
(nuclear-scale) distances. Specifically, the effective coupling constant of perturbative 
QCD increases at low energies, becoming of order 1 at energies ~ 200 MeV (conversely, 
the coupling constant approaches at short distances or high energies, a property called 
"asymptotic freedom"). Therefore the main analytical tool used to study quantum field 
theory - perturbation theory - breaks down as this energy scale is approached from 
above. 

Qualitatively, we expect QCD below this energy scale to "confine", or tightly bind 
quarks into color-neutral bound states, which are the familiar hadrons of particle physics 
(such as protons, neutrons, pions and other particles). The analogous theory without 
quarks (non-Abelian gauge theory, also called Yang-Mills theory) is also asymptotically 
free and is expected to manifest similar behavior at low energies: the massless gluons of 
Yang-Mills perturbation theory, which mediate the strong nuclear force, bind together 
into hadronic "glueball" bound states and become massive. Approximate numerical 
results in QCD and Yang-Mills theory (such as the value of the hadron masses) can be 
obtained by simulating the theory on a discrete spacetime lattice, and various qualitative 
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proposals have been made for the mechanism of confinement, but a solid theoretical 
understanding of confinement is still lacking 1 . 

In the absence of analytical tools for studying non-perturbative phenomena in QCD 
such as confinement, one alternative is to turn to related models in the hope of finding 
a more tractable problem that may nonetheless provide insight into the theory of inter- 
est. A profitable tool in this regard is supersymmetry, a symmetry that relates bosonic 
and fermionic degrees of freedom. The extra symmetry constraints present in the super- 
symmetric version of Yang-Mills theory and QCD are surprisingly tight and allow for 
greater depth of analytic computation; at the same time, the super symmetric versions 
of Yang-Mills and QCD are expected to share many of the same qualitative features, in 
particular confinement at low energies. 

In fact, there are several theoretical and experimental indications - and widespread 
anticipation among high energy particle theorists - that supersymmetry may be realized 
in nature at suitably high energies. Thus, the study of supersymmetric gauge theories 
may be directly relevant for describing the nature of fundamental interactions at suffi- 
ciently high energies. 

It has long been suspected that four-dimensional gauge theories such as QCD are 
related to string theories. The tube of confined gauge field flux that extends between 
two quarks has string-like properties, and in fact, modern string theory emerged from an 
attempt to model the strong interactions. However, despite over three decades of inten- 
sive study there is still no known consistent quantum theory of strings propagating in 
four dimensions; for example, worldsheet anomaly cancellation of the supersymmetric 
string requires the (suitably generalized) dimension of spacetime to be 10. 

The resolution to this dichotomy is that four-dimensional gauge theories may be 
equivalent to (limits of) string theories in higher dimensions; the dynamics of strings 

'See | ApoAD| for a classical problem of comparable difficulty. 
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propagating in the extra dimensions can give rise to gauge dynamics in four dimensions. 
These "gauge/string dualities" have provided many fascinating and unexpected results, 
some of which are the subject of this thesis. 

As in the heuristic example of QCD, open strings carry matter degrees of freedom 
("quarks") at their endpoints, and their excitation spectrum contains a massless spin-1 
particle. Thus, open strings give rise to matter coupled to gauge fields. Taking into 
account string interactions, the endpoints of open strings may join together to form a 
closed string. Closed strings include a massless spin-2 particle in their excitation spec- 
trum; this particle must couple to the stress-energy tensor of the theory, and the space- 
time theory is required by consistency to have diffeomorphism invariance. The spin-2 
particle is therefore identified with the graviton, and quantum theories including closed 
strings are theories of quantum gravity that reduce in the classical limit to classical gen- 
eral relativity coupled to additional fields. Thus, string theory has the potential to unify 
the interactions of matter with all four fundamental forces in a consistent quantum the- 
ory; this is a long-standing theoretical problem that has resisted many previous attempts 
at solution. 

The link between these two aspects of string theory, and the main string-theoretical 
tool for studying gauge theories in the modern context, are D-branes. These are extended 
"membrane" objects, of various dimensions, which are required by non-perturbative 
consistency to be present in the spectrum: when the theory contains open strings, these 
strings may only end on a D-brane. Therefore, the matter fields at the endpoints of the 
strings are confined to live on the D-brane, and open strings with both endpoints on 
the brane give rise to gauge fields propagating along it. Thus, the study of D-branes 
and strings propagating in appropriate 10 dimensional geometries can teach us about 
four-dimensional gauge theories. 
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This is one of the great advantages of string theory, fully realized only since the mid 
1990s; it can be used to translate certain problems of quantum field theory into geomet- 
rical language. This allows the application of powerful geometrical tools to study the 
corresponding quantum field theories. In cases where theoretical control of the calcula- 
tion is presently available, the corresponding field theories are typically supersymmet- 
ric, and the more supersymmetries that are present, the greater the constraints on the 
mathematical structures that underly the theory. 

String theory is now known to possess many remarkable properties, and while there 
remain many difficult problems to solve before it can be quantitatively applied to study 
the physics of our observed universe, it has nonetheless provided deep insights into 
many aspects of theoretical physics and mathematics. In this thesis, we will describe a 
set of tools that have emerged from string theory over the past few years, which allow 
the computation of exact results in a class of confining supersymmetric gauge theories 
at low energies. These string theoretical tools have provided some unexpected insights 
into the structure of quantum field theory. 

Central to the analysis is a particular class of hyperelliptic curves related to a string 
theory background geometry, the periods of which encode the superpotential of the 
gauge theory and define its vacuum structure. These "spectral curves" also emerge from 
the study of a number of mathematical systems that appear at first sight to be unrelated 
to the gauge theory (such as matrix integrals, and integrable systems), and understand- 
ing this connection provides new insights into the structure of the vacua of the quantum 
field theory 2 . 

2 Conversely, this relationship provided a link between previously unrelated areas of mathematics, 
for example that the combinatorics of planar diagrams is related to the special geometry of Calabi-Yau 
manifolds. 
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To provide context for the later results on supersymmetric gauge theories, we will 
begin by reviewing some known techniques and results on the vacuum structure of non- 
supersymmetric gauge theories. We will explain the limitations of these calculations, 
and describe how they are avoided in supersymmetric theories. The remainder of the 
thesis will discuss various techniques that have emerged from string theory and allow 
the computation of exact results about the low energy structure of supersymmetric gauge 
theories. 

This thesis is based on material previously published in the original collaborative 
works «ACH+03llKW03i and on the review article [Ken04|, although some details and 
aspects of the composition are new. 
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Chapter 2 

Effective Potentials in Quantum Field 
Theories 

When a quantum field theory possesses continuous symmetries, the form of the effec- 
tive potential (the non-derivative terms in the effective Lagrangian) is constrained by the 
corresponding (anomalous) Ward identities, which give rise to partial differential equa- 
tions that must be satisfied by the quantum corrected effective potential. For example, as 
we will discuss in section 12.1.21 the differential equation associated to the anomalously 
broken scaling symmetry is the Callan-Symanzik equation. 

The background field method can be used to derive the one-loop effective action from 
the path integral of the theory; in theories with non-trivial vacua, such as asymptotically 
free theories, this gives an approximation to the vacuum state. Evaluating the 1-loop 
effective action is equivalent to the summation of an infinite class of Feynman diagrams 
where one includes the couplings of a set of fluctuating fields to a classical background 
field, but ignore the self-interactions of the fluctuating fields. 

We begin by studying the Gross-Neveu model, a two-dimensional theory of chiral 
fermions which is asymptotically free. This model exhibits several of the features of 
more interesting four-dimensional theories such as Yang-Mills theory and QCD, includ- 
ing asymptotic freedom and spontaneous chiral symmetry breaking. We will solve for 
the 1-loop effective potential of this model, as a warm-up exercise for studying four- 
dimensional gauge theories. 
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Due to the Landau pole (divergence of the perturbative gauge coupling at low ener- 
gies), the one-loop approximation to the Yang-Mills effective potential cannot be extrap- 
olated to the vacuum of the theory, but it gives a qualitative picture of some of the 
features of the vacuum. When the theory has supersymmetry, the constraints on the 
effective (super)potential become much more powerful, and the one-loop perturbative 
gauge theory computations can be extrapolated all the way to low energies to obtain 
exact, non-perturbative information about the vacuum. 



2.1 A toy model: the Gross-Neveu model 

The Gross-Neveu model IGN74II is a simple model that exhibits spontaneous symmetry 
breaking through a quantum-mechanical symmetry-violation. It is a two-dimensional, 
asymptotically-free theory of N massless interacting fermions, with Lagrangian: 

£ GiY =?^ + y(^) 2 (2.1) 
The classical Lagrangian has a discrete chiral symmetry 

A -> l5ipi ipi -> -V'iTs (2.2) 

By summing the contribution of Feynman diagrams with vanishing external momenta, 
we will derive the effective potential of the Gross-Neveu model, and find that the chiral 
symmetry is spontaneously broken in the quantum theory. This perturbative 1-loop 
computation provides exact non-perturbative results about the vacuum of the theory at 
large N. 

A useful technique for studying the response of quantum field theories to non-trivial 
field backgrounds is the background field method. One splits the external field into a 
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classical, background field, and a fluctuating quantum field, and then evaluates the path 
integral perturbatively in the fluctuations around the given background. We will make 
full use of the background field method when we study non-Abelian Yang-Mills theory 
in section 17.2 .21 Because this technique is non-perturbative in the background field, it 
can be used to probe for phenomena that are invisible in perturbation theory around the 
usual zero-field background. 

Fermionic (Grassman- valued) fields are not usually considered as classical field the- 
ories, for example as possible background fields for a quantum field theory calculation. 
However, fermionic quantum fields can pair up and form a composite bosonic field 
a ~ ipip which can attain a vacuum expectation value. The Gross-Neveu Lagrangian 
can be rewritten as 

t = ipitfiipi - ^cr 2 - aip^ (2.3) 

which re-expresses it in terms of a coupling to the composite bosonic operator a. This 
field is treated as a non-dynamical, external background field since it has no kinetic 
term. It is easily verified that integrating over this auxiliary o field recovers the original 
form of the Lagrangian (12. 1 b . The Feynman rules for (12.31) are shown in figure |2~T1 

We will analyze this theory in two ways: by performing a path integral computation 
that amounts to summing the Feynman diagrams that can contribute to the effective 
potential of the theory due to the interaction with the external a field, and by using the 
anomalously broken scale invariance to constrain the form of quantum corrections to the 
potential. 
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Figure 2.1: The Feynman rules for the Gross-Neveu Lagrangian (12.31) 

2.1.1 Path-integral computation of the effective potential 

We can probe the response of the Gross-Neveu model to the formation of a non-zero 
fermion condensate by introducing an external source J for the field a into the path 
integral, finding the minimum-energy field configurations in the presence of the source, 
and then turning off the source. Define 



where — E[J] is the generating functional of connected correlation functions of a. 
Define the classical field 



the vacuum expectation value of a(x) in the presence of the source J. Then the Legendre 
transform of the energy functional — E[J] defines the effective action r(<r c /) 



e -*E[J] 



Va J] J] Vij) i exp (z (C(a, ^, ^) + J a) ) (2.4) 



5E 
J7 



(Q\a(x)\0)j 



(2.5) 




(2.6) 
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subject to the constraint 

ST(a d ) 



-J(x) (2.7) 



5a d (x) 

Thus, turning off the source J we obtain that the stable configurations for the external 
field a d are those for which 

T¥i = ° < 2 - 8 > 

da d {x) 

In the translation-invariant vacuum states of the theory, a d (x) is constant, and the effec- 
tive action can be written as 

r[a d ] = -(VT)V ef (a d ) (2.9) 

where V is the 3-dimensional volume, T is the time interval of the integration region, 
and we defined V e fda d ) the effective potential for the classical field o d . The vacua of 
the theory satisfy 

dVff{a c i) 

-f = (2.10) 

aa d 

The effective action is the generating functional of 1 -particle irreducible (1PI) cor- 
relation functions of the a field. Therefore in the background of a d 

V eS (a d ) = ^ F "(°' °> • • • > °) C 2 - 11 ) 

where the 1PI diagrams that contribute to T n carry external momenta on all legs, and 
each leg comes with a coupling to the background field. To 1-loop order, the diagrams 
contributing to the effective potential are shown in figure 12.21 Since they all involve 
a single fermion loop, we can evaluate the 1PI generating functional to 1-loop order 
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o • o • O 



Figure 2.2: Diagrams contributing to the 1-loop effective potential for the background 
field a 

by integrating over the fermions, which appear quadratically in the path integral of the 
original Lagrangian: 



Vae Sd2x ^ det(i0 + a) N 



J Vae^ d2xC{(7) (2.12) 



with 



a 2 



C(a) = -— + tN\ogdet^ + a) (2.13) 

Using the two-dimensional gamma matrices 7 = a 2 , 7 1 = 1a 1 and performing a 
Fourier transform, we can evaluate the determinant in (12.131 ): 



log det(ift + a) 



d 2 p 



log det(iA + cr) 



(27T) 2 

d 2 p 
(2ttV 

\ tp + ipi 



log det 



a 



-ip + ipt 
a 



d 2 p 



log(cr - p ) 



(2.14) 
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Therefore 



= + *N I 7^ log(a 2 - p 2 ) (2.15) 



2g* J (2tt) 2 

The 1-loop 1PI correlation functions of (12.31 ) may be derived from £(cr), and to this 
order we can identify the Lagrangian C(a) with the effective Lagrangian associated to 
the effective action T(a) = J d d x C £ jj{a), or in other words 

%M = -C(a d ) (2.16) 
We can recover the diagram sum explicitly by writing 

log(a 2 -p 2 ) = log(-p 2 (l-^)) 

= log(l - ^f) + log(-p 2 ) 

pZ 



logfl 



M 2 . 

p 2 ■ 



oo ^ 

1 <™^2n 



E-(-) 

^ n p 



n=l 

oo 

TrVi( , 

n=l 



ta)i) 2n (2.17) 



where in the third line we dropped the second term since it just gives rise to an infinite 
constant upon Wick rotation and integration over p. Comparing to the Feynman rules in 
figure l2~Tl each term in the series corresponds to a 1-loop diagram of the form shown 
in figure l2~2l therefore, integrating over the fermions to quadratic order is equivalent to 
computing the 1-loop diagram sum to all orders. 

Returning to the 1-loop effective Lagrangian, the integral (12.151) is divergent and 
needs to be regularized. Wick rotating to Euclidean space and using dimensional regu- 
larization we obtain 
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2,f^ J (2ir)"da VpI + W 

- ¥ + "iiSFS^ r(-a) UJ JJ o 2 ' l8) 

where A = a 2 . Using the expansion of T{x) near its poles, T(x) ~ J(" x + n ) — 7 + 1 + 
...+ - + 0(a; + n) and T(x + 1) = xr(a;) we expand T(— a — |) and write the singular 
terms in the form suitable for the modified minimal subtraction scheme (adapted to 2 
dimensions): 



r(i-f) f-i\ l - d/2 1 (1 



(4^/ 2 \Aj ^ ^ Ve -7 + log4.-logA + 0(e ) 



1 1 f A 

47T yU 



log(-) (2.19) 



We obtain for the effective potential 



log-f-l (2.20) 



4tt V & A 2 



where in the second line we defined the dynamical scale A 2 = /i 2 exp(-^p-). The poten- 
tial (12.201) is of Coleman-Weinberg type ICW73I and has the form shown in figure 12.31 
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Extremizing (I2.20I ). we find that what was the classical minimum (a d ) = is now 
a local maximum, and there are degenerate global minima at (ad) = ±A. Thus, the 
original "perturbative" vacuum can minimize its energy by spontaneously generating a 
background of paired fermions, 

<^> = ^> = ±4 ex p(^) ( 2 - 21 ) 

9 9 Ng 2 

and since this fermion bilinear does not respect the chiral symmetry (12.21) . the Gross- 
Neveu model exhibits spontaneous chiral symmetry breaking. 

Higher loop corrections to the effective potential necessarily involve a propagators 
and are therefore suppressed by powers of g; in fact all higher loop corrections van- 
ish in the 't Hooft limit iV -> oo, g -> 0, g 2 N = const. IGN74I . Therefore, in this 
limit the 1-loop result is exact. Unfortunately, for most interesting non- super symmetric 
theories (such as Yang-Mills or QCD) the higher-loop corrections do not vanish in this 
limit, and the infinite diagram series cannot be summed explicitly even at large N l ; 

'Although a generating function that enumerates the infinite series of Feynman diagrams is known for 
QCD |tH99l 
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the complication comes from performing the loop momentum integrals at higher orders. 
However we can obtain partial results by organizing the diagrams as a loop expansion: 
in section IT. 2 .21 we will show how summing the one-loop diagrams for Yang-Mills the- 
ory in a covariantly constant background field strength gives a (not particularly good) 
approximation to the vacuum state of Yang-Mills theory. 

However, simplifications even more powerful than those of the Gross-Neveu model 
were observed recently in certain four-dimensional J\f = 1 theories, where supersym- 
metry provides additional constraints on the effective potential that allows us to sum the 
diagram expansion to all orders. We will come back to this in section l2~4l 

The value of the fermion condensate (ipip) = ±p-exp(-^) is a non-perturbative 
quantity, since its Taylor expansion around g = vanishes to all orders. Therefore the 
non-trivial vacuum of the Gross-Neveu model is invisible in the perturbation theory of 
the original Lagrangian (12.11) . which preserves chiral symmetry to all orders. It was only 
by rewriting the Lagrangian by introducing a coupling to the appropriate background 
field that we could probe for the existence of a chiral symmetry breaking condensate. 
We have seen that by introducing an appropriate variable in which to perform a perturba- 
tive loop expansion (the composite background field a), we can obtain non-perturbative 
information about the vacuum of the theory, order by order in the perturbative evaluation 
of a different Lagrangian. 

2.1.2 Anomalous symmetries and effective potentials 

In quantum field theories, continuous symmetries of the classical Lagrangian may some- 
times be violated in the quantum theory. An example of an anomalous symmetry are 
scale transformations (dilatations) in massless field theories 2 . The continuous dilatation 

2 Another anomalous symmetry is axial rotations of massless Dirac fermions in gauge theories; the cor- 
responding effective Lagrangian including quantum corrections from the axial anomaly can be obtained 
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symmetry is associated to a current = Q^ v x u , where is the stress-energy tensor 
of the theory, defined by 



Classically the dilatation current is conserved; d^D^ = 6^ = 0. However under a 
change of renormalization scale this symmetry is broken by the running of the coupling 
constant (see KPSIO . and the one-loop trace anomaly is given by: 



The trace anomaly receives contributions from all orders in perturbation theory, as well 
as possible non-perturbative corrections, through the beta function. 

In a quantum field theory the "charge" of fields under a scale transformation (their 
scaling dimension) may receive quantum corrections as we change the renormaliza- 
tion scale; operators can have anomalous dimensions. The Callan-Symanzik equation 
encodes the scaling behavior of the effective potential under a change of renormalization 
scale (renormalization group invariance): 



where d is the space-time dimension, dj are the classical scaling dimensions of the oper- 
ators Oi, are their anomalous dimensions, and /i is the renormalization scale. This 
equation imposes that the effective potential must scale with dimension d, and reproduce 
the trace anomaly under a scale transformation. 

by similar techniques, and has been used to study the role of the anomaly in the low-energy dynamics of 
mesons IDVV80IIWit80l . 




(2.22) 



dpD* = P{g)—C 



(2.23) 



d-J2(d i + ^ i )o i — + m - + ,- v eff =o 



(2.24) 
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In order to use the Callan-Symanzik to obtain predictions about the form of V e m 
we need to know the (3 function and anomalous dimensions 7. These are typically only 
known through explicit loop calculations, such as the one we did in the previous section. 
However, as we will discuss in section |2~.2.31 once we know j3 and 7 from a particular 
calculation, we can use the Callan-Symanzik equation to constrain the allowed form of 
the effective potential for an arbitrary field background. 

We impose 



2 - (1 + 7 >l + ^ ) | + 4j 



V eff =0 (2.25) 



and find that (3(g) = — 4^-, 7 CT = 3 . As we noted in the previous section, in the 't 
Hooft limit these quantities are exact. 



2.2 Four-dimensional gauge theories 

Before considering non-Abelian Yang-Mills theory, it is instructive to review the calcu- 
lation of the effective potential for QED in external electromagnetic fields, which shares 
many technical features with the Yang-Mills case. These results were first obtained by 
Euler and Heisenberg in 1936 [HE36|, and were cast in a rigorous quantum field the- 
ory framework by Schwinger in 1951 HSch51l . The presentation here includes elements 
from llSS75l|Fly80| . 

3 The field a has vanishing anomalous dimension due to the normalization of the Lagrangian (12.31 . A 
wavefunction renormalization of a cannot be balanced by a coupling-constant renormalization since the 
coefficient of the a interaction term is fixed to 1. 
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2.2.1 QED 

The Lagrangian of QED is 



£ = —^vF^ + ijpij + mW (2.26) 

where the covariant derivative is = <9 M + ieA^. As in the previous section, the 
effective action for the gauge field is given to 1-loop order by 



= &et(ip-m)e-^ dAxF2 

= exp(i J d A xC e ^) (2.27) 



where we defined the 1-loop effective Lagrangian 



% = —F^F"" -ilogdet(ip-e/L-m) 

(2.28) 

For comparison to Yang-Mills theory in the next section, we henceforth restrict to 
massless electrons, although the massive case can be easily treated in a similar manner. 
To evaluate the fermion determinant det (ip) it is convenient to evaluate the determinant 
of {ip) 2 and take the square root. Expanding and using the anticommutation relation 
= we find 



= —D 2 - C -a, u F^ 



(2.29) 
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where |[7 M , 7"] 



cr^ is the generator of Lorentz transformations on the spin-^ repre- 



sentation. Therefore 




e 



r ) 



(2.30) 



As we discussed in the previous section, the determinant corresponds to summing 
up the infinite series of 1-loop Feynman diagrams of the theory, where the electron runs 
in the loop, and we consider arbitrary insertions of the background gauge field. The 
one-loop effective Lagrangian for massless QED is therefore 



This Lagrangian exhibits the anomalous magnetic moment interaction ^a^F^ of the 
electron with the background electromagnetic field. A similar magnetic moment inter- 
action for the charged gluons of Yang-Mills theory will be vital for understanding the 
vacuum properties of that theory. 

In diagonalizing this operator one needs the eigenvalues of the field strengths F^. 
Defining the Lorentz scalar and pseudo-scalars 



%= -\F^ U ~ ^Trlog((p, - A,f - e -a, v Fn 



(2.31) 



T 



-F^ = -{&-E*) 




(2.32) 



where F^ 



^ie^ upa F pa is the dual field- strength tensor. Using the identities 
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F^F pv = -5;g (2.33) 

F W FP V - F w Ff> v = 25;? (2.34) 

the eigenvalues A of are found to satisfy 

A 4 + 2T\ 2 -G 2 = (2.35) 
which has solution ±A^, ±A^, with 

A« = -L ((jF + ^)V2 + ( JF _^)l/2 ) (236) 

A( 2 ) = _L((jr + ,g)i/a_(jr_^)i/2) (237) 
v 2 

The magnetic moment operator satisfies 

(i^FH 2 = 2(^ + 75 £) (2.38) 
therefore using 7! = — 1 and (12.361 ) the eigenvalues are 

±{2{F ±iQ)Y' 2 (2.39) 

In a particular Lorentz frame, a constant magnetic field may be specified by taking 
Q = 0, T > 0, and the eigenvalues A are real. For a constant electric field Q = 0, T < 
they are purely imaginary; this difference is the cause of the vacuum instability we will 
find for the constant electric field. 
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First, we consider a constant magnetic field, which we take to be along the 3 direc- 
tion, A = (0, 0, -Bx u 0), B > 0, and we have Q = 0, T = \B 2 , and 

\<r^ v ~ diag( eJ B, eB, -eB, -eB) (2.40) 
In this gauge the d'Alembertian D 2 becomes 

D 2 = pi - v\ ~ (P2 + e^i) 2 - p\ (2.41) 
and becomes after a unitary transformation 

£)2 = eVlPa /eB^ -p\- e 2 B 2 x\ - p 2) e -*PiP2/eB ( 2 .42) 

where we have used the commutation relations [a^, p„] = ig^ v , and in particular 
[ Xl ,e api ] = iae api . 

Therefore the 1-loop contribution to the effective Lagrangian is 

£i = -iTrlog (e^ eB (pl - V \- e 2 B 2 x\ - pl)e^' eB - ^F^) (2.43) 

To evaluate this trace, we use the identity 

-i e f°° 

log(x) = lim — dtt- 1+e e~ ltx (2.44) 

e^o r(l + e) J 

This is related to the method used by Schwinger [Sch51 1 (who introduced a lower cut- 
off into the integral instead of dimensionally continuing the argument), and amounts 
to rewriting the four-dimensional space-time loop momentum integral as the world-line 
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integral of a particle moving in an external potential. This is a close analogy of the 
world- sheet formalism of string theory; the world-line proper time parameter t corre- 
sponds to a "world-line modulus" of the loop in the Feynman graph 4 . 




oo 



2r(l + e) 



Tr 



dt t~ l + e e W2Pi/eB e -it(;pl-pl 2 -e 2 B 2 xl-pl) e -i V2 pi/eB e it^a^F^ 



2T(l + e) 




A=±l 

For suitably large e the integral converges, therefore this representation regulates the 
calculation. In the second line we evaluated the trace over the anomalous magnetic 
moment operator using (I2.40I) . since the operator commutes with everything else in the 



The an alogy between string theory and the Schwinger formulation of loop integrals was used in 
|DGL + Q3l to calculate effective superpotentials in theories with Af = 1 supersymmetry, by reducing 
a topological string theory calculation to a field theory calculation in Schwinger's formalism. We will 
explain some key points of this work in section lZ4l 
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trace. The remaining trace may be evaluated as follows iflzll 

d A pd A V ' (xl^^le^/^e-^-P^ 62 ^^^^^-^ 1 / 65 !^)^'^) 

d 4 p d 4 p' 6l< * P ^ (p| e 1 P2Pi/ eg e -» f (pg-P?- e2fl23! ?-^) e -^Pi/«g|^) 
(2tt) 4 

(fp d*p' P — c ^(P2Pi)/ej?-(p^P , 1 )/eB) e -rt(pg-pj) 

(2tt) 

(pi|e- , -*W- afla ^bi>(5 s ((p-j/)o, W ) 

^ e - l t( P g-pI) el ^-o;')(^+P2/eB)^| el t(pf+e^^?)| a; /^ 



(2vr) 3 



(2 7 r) 2 (zt) 1 /2(_ z t)i/2 



da; da;' % - a/)(u;|e^I +e B 



exp(z£(n + -)2eB) (2.46) 



(2vr) 2 £^ ^ v v 2' 

v ' n=0 

where we used the result for the energy levels of a harmonic oscillator 



Trexp(z£(^- + ^Q 2 )) = f^exp(tt(n + (2.47) 

n=0 

Therefore the effective Lagrangian reduces to 



£l = ^wT^T / rftr2+e S exp(ze££A) V exp(z eJ B£(2n + 1)) 
8vr 2 r(l + e) y ^ 



A=±l n=0 



{2eB) 1 - £ dtr 2+e -—^ (2.48) 



87T 2 r(l + e) v ' 7 1-e 
Rotating the integration contour £ — > z£ we obtain 



£ x = r(2e5)- £ / d££" 2 + e ^^ (2-49) 

47r 2 r(l + e) V Wo 
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The integral may now be evaluated using the identity 



e -vt 



dtt- 1 — - = T(a)C(a,u) (2.50) 

1 — e 1 

where ((a, v) = Yjn=o{ u + n )~ a ls tne generalized Riemann zeta function. Therefore 



e 2 B 2 ( 1 V-IY-I 



(2.51) 



In taking the limit e — > 0, we renormalize the expression using a variant of the MS 
scheme [PS l 5 : 



m /1V 1 log f 4' I (2.52, 



(47r) 2 + e V2A7 4tt 2 ° \/j, 2 

and use the property of the ^-function 



C(-m, i/) = — — (2.53) 

m + 1 



where m = 0, 1, . . ., and B m+ i(u) are the Bernoulli polynomials, in particular 
E>2(x) = x 2 — x + 1/6. Putting this all together, we find for the effective potential 



V eJf =^- e ^log(eB/v 2 ) = ^-^log(e5//i 2 ) (2.54) 

3 

where we recognize the 1-loop QED /3-function coefficient b = ^l-?- This potential is 
plotted in figure l2~4l 



5 The difference is that we also subtract the log 2 coming from the coefficient of A 
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Figure 2.4: 1-loop effective potential for QED in a constant background magnetic 
field. The apparent instability at large magnetic field strengths is an artifact of the 
1-loop approximation. 

For small external fields eB < /j 2 the second term is positive, and the effective 
potential has a local minimum at B = 0. At larger field strengths there appears to be a 
local maximum and the potential eventually becomes arbitrarily negative. However, in 
precisely this limit the 1-loop approximation breaks down, because the quantum correc- 
tion term dominates and is no longer small compared to the classical term. Therefore, 
for large enough magnetic fields one needs to also consider the higher loop corrections. 

We turn now to the electric case. Using the form of T = \[B 2 — E 2 ), we may obtain 
the effective potential for a constant background electric field E ^ 0, B = by formally 
continuing B — > iB = E in (12.541) . This introduces a factor of % into the argument of 
the logarithm, and therefore the effective Lagrangian in a background electric field is 
complex. 

Since the amplitude for a vacuum in the far past to remain in the far future is given 

by 

(0 + |0_) = e ir (2.55) 
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the probability of vacuum decay, per unit time and volume, is given by 



2Im£ 



Im 



127T 2 



log(i) 



24tt 



(2.56) 



and the constant electric field background is unstable against pair production of 
positron/electron pairs. 

The result for a non-zero electron mass can also be computed following the above 
steps, and one finds 



which is non-perturbative in the RG-invariant field combination eE. Again we see that 
the background field method produces non-perturbative information from a perturbative 
calculation. 

In a general constant background with T 7^ 0, Q ^ 0, the effective Lagrangian is 
that of Euler and Heisenberg IH E36L which takes the form (before regularization and 
renormalization) 

-C 1 = A r dt / e2a ^osh(eat)cos(e6t) \ 
871"^ J \ sinh(eat) sm(eot) / 

where a 2 — b 2 = E 2 — B 2 , ab = E ■ B. A list of references to recent work on this 
Lagrangian and related matters may be found in [Dun04|. 

2.2.2 Yang-Mills theory 

To calculate the 1-loop effective action for four-dimensional Yang-Mills theory we again 
use the background field method. This calculation and related results were developed by 




(2.57) 



N=l 
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a number of authors, including IDR M751 IBMS77L IMS781 IN0781 IPT781 IBW791 lYC?80l 
FIygUlFPgTllJWZ81ll . 

The Yang-Mills Lagrangian is 



£ = -\F; v F^ a (2.59) 

We split the gauge field into a classical background field A and a fluctuating quantum 
field a: 



Afr) -> A;(x) + a;(x) (2.60) 

The covariant derivative (D^) ac = d^5 ac + igf abc A b ^ is defined with respect to the back- 
ground gauge field, and we will integrate over the quantum field a in the path integral. 
Then the field strength becomes 

F% - F % + D »< ~ D K + i9f abc ay v (2.61) 
In background gauge D^A^ a = 0, the gauge-fixed Lagrangian is 



£ = -J(f; + m-M + ^A;<) 2 

--(5V) 2 + c a (—(D 2 ) ac - igD» ad f dbc al)c c (2.62) 

where c, c are the Faddeev-Popov ghosts corresponding to the gauge fixing. 

As before, the effective action to 1-loop order is given by evaluating the path integral 

<™ = JV**D#1<« (2.63) 
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to quadratic order in the fluctuations. Expanding (12.621) to quadratic order, we find 



tquad = ~\< l(-D 2 ) a V v ~ 2igF^ c r b ] fl J + c a [—(D 2 ) ab ] c b (2.64) 



As in QED, the new interaction term — 2igF ,luc f co6 o /io a is an anomalous magnetic 
moment interaction of two spin-1 gluons with the background field F^ c . Introducing 
the generator of spin- 1 Lorentz transformations 



the operator —2igF puc f cab can be rewritten as —2i(^F^ a J pcT ) fMU f cab , emphasizing the 
similarity to the operator (12.291) for spin-| electrons in QED. The spin interaction for 
the ghost fields vanishes since they have spin 0. 

Therefore the path integral to 1-loop order is Gaussian and can be evaluated, giving 
the 1-loop effective Lagrangian 



We can evaluate these determinants by restricting to covariantly constant fluctuations of 
the gauge fields: 



(J pa U = «0M - «) 



(2.65) 




DpF^ = <£> [D p , F^} = 



(2.67) 
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where we write the field strength as a matrix in colour space (F pu ) ab = f abc F^ u , and the 
second form follows because the covariant derivative in the adjoint representation acts 
by matrix commutation. Using the Jacobi identity [D a , [D p , F pu ]] +perm. = it follows 
that 

[F pv ,F pa ]=0 (2.68) 

i.e. the colour matrices F pv form a commuting set and may be simultaneously diago- 
nalized. In other words, by a gauge transformation we may rotate a given gauge field 
configuration into the Cartan subalgebra. Then 



£} = -TrJ2 l °s(-D {a)2 g^ -2igF^) -j^Trlog((-D (tt) ) 2 ) 

a a. 

= jTr^logt-D^ 2 ^-^^) -2z^Trlog((-.D (a) ) 2 ) (2.69) 

where the sum is over the positive roots a of Q. In the second line we used that each 
root a is paired with a negative root —a, and the zero roots do not contribute. We also 
defined effective quantities 



£>W = dp + igajAi 

P( a ) — (y F j 

Af = a,-4 (2.70) 



in terms of the simple roots (a 1: . . . , ct r ), r = rank(^), which span root space. 

In other words, we have reduced the computation of the 1-loop effective action for 
a non-abelian gauge group Q to that of an Abelian U(l) r gauge theory, where the j'th 
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"photon" carries charges gctj with respect to the different £7(1) gauge factors. The situ- 
ation is therefore quite similar to that of QED, which we studied in the previous section, 
except there is more than one type of "electromagnetic field", and the charged particles 
are spin-1 photons, not spin-| electrons. 

At this point we need to choose the orientation for the effective £7(1) gauge fields 
in four-dimensional space; when the rank of the gauge group is larger than 1, the "elec- 
tromagnetic fields" may point in different spatial directions. Most of the early work on 
this problem either considered SU(2) IB MS771IMS 78I. or chose to align all the effective 
£7(1) gauge fields parallel 1AN07 9I. However, it was shown in subsequent work that 
for 2 < N < 4 the lowest-energy configuration is to choose the fields to be mutually 
orthogonal [Fly 80 1 . ForiV > 4, i.e. rank higher than 3, it is no longer possible to choose 
all vectors to be orthogonal in three-dimensional space, and for N — > oo the minimum 
energy configuration corresponds to an isotropic distribution in space IF P8T1IJWZ81I . 

For simplicity, we will henceforth restrict to the SU(2) case. The essential features 
are seen in this case; in particular we will see that any choice of covariantly constant 
field strength gives rise to a vacuum instability, and therefore the 1-loop result is at best 
only an approximation to the true vacuum. This instability persists for the non-parallel 
gauge field orientations mentioned above. 

We can now proceed as in section 12.2.11 Again taking a constant magnetic field, 
the eigenvalues of —2%gF^ v are (±2gB, 0, 0). The two zero eigenvalues cancel with the 
contribution from the ghost determinant in (12.69b . giving 

C 1 = Trlog(-/J 2 - 2XB) (2.71) 
x=±i 

After manipulations similar to QED, we find 
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C 1 = J dtr^x £ exp(-it(l - 2A)) 



A=±l 



r(i + e) y 

^exp(-2z£A0 (2.72) 



oo 

X 

N=0 



Note that we can no longer unconditionally rotate the contour by taking £ — ► — z£, 
because the mode with (A, AT) = (1,0) would diverge like e*. This is the unstable 
mode found by Nielsen and Olesen HN078L which will give rise to an imaginary part 
for the 1-loop effective Lagrangian even in the magnetic case. To proceed, we subtract 
and add the (A, N) = (1, 0) term: 



3tt 2 ^ ' r(l + e) i J l-e~ l/ 

POO 

+i 1+e / dtt~ 2+e e l 



e 



-{gB)- £ — <^ i 1+e / d££ 



8vr 2 ^ y r(l + e) [Jo 1 - e" 

/■oo 

/ dtt~ 2+t e l 



o 



8tt 2 ^ y T(l + e) 1 V 2 / 7n 1-e"* 

/•oo 



8tt 2 ^ y r(l + e) i V 2 / 7o 1-e"* 

+(-l) e d*r 2+e e-*| (2.73) 

where we rotated the two integration contours by £ — > — z£, £ — > z£ respectively. The 
integrals may now be evaluated in terms of zeta functions, giving 
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1 




*<- i+e '!» + S) 




2<? 



M^/V ) -2(C(-l, \) + C(-i, §)) + i 



(2.74) 



where in the last line we recognized the 1-loop /3-function coefficient. As before the 
pure electric field result may be obtained by analytic continuation. If we consider a 
background with Q ^ 0, then the effective Lagrangian will be a generalization of the 
Euler-Heisenberg Lagrangian (12.581) [BMS77, MS78|. In all cases the background is 
unstable, in contrast to QED, for which only the electric background is unstable. 

Note that because of asymptotic freedom the sign of the 1-loop term is opposite to 
that of QED (12.541) : therefore the effective potential has a similar form to figure 12.31 
The lesson we can draw from this analysis is that the "perturbative vacuum", where we 
consider excitations around the zero-field background, is an unstable field configura- 
tion. The Yang-Mills vacuum lowers its energy by spontaneously generating a non-zero 
background field. This can be seen as a vacuum anti-screening effect by the gluons, 
which are charged under the gauge group and can act as sources for other gluons. Turn- 
ing on a covariantly constant background field indeed lowers the vacuum energy, but 
this field configuration is itself unstable (not to mention violating Lorentz invariance), 
so the "true" vacuum is some other field configuration. An ansatz for the vacuum (the 
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"Copenhagen vacuum") was proposed in HNN79L based on exciting the unstable mode 
of the constant-field vacuum. 

The background field method is non-perturbative in the background field (since it 
is not used as an expansion parameter), which allowed us to make some progress, but 
excitations around this field still must be calculated perturbatively. This means that 
we can only trust our 1-loop calculation when the effective coupling constant is small, 
however this is counteracted by the negative sign of the 1-loop /^-function, which tells 
us that g will grow towards the IR. 

Explicitly, to 1-loop order the running of the Yang-Mills coupling constant is given 

by 



9 2 ef f{q) = u f - - - (2-75) 

37 l + ggfelog(9/^) 



which diverges at the finite energy scale 



96vr 2 iV\ 

n? 



g = /xexp( 7t^-) = Aym ( 2 - 76 ) 



Therefore, we can not trust our 1-loop effective potential at energies comparable to or 
lower than AyM- Nevertheless, it is expected (based on lattice simulations and other the- 
oretical work) that the qualitative picture remains true, and the vacuum of Yang-Mills 
theory is associated to non-trivial gauge field backgrounds, which give rise to confine- 
ment, generation of a mass gap (the appearance of massive glueballs in the spectrum 
replacing the massless gluons), and other poorly-understood low-energy physics. 

2.2.3 Constraints on the effective potential from the trace anomaly 

We have seen that the effective potential of quantum field theories must be consistent 
with the trace anomaly, in particular it satisfies the Callan-Symanzik equation. Once 
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we have calculated the quantities (3 and 7 for a particular theory, we can use the Callan- 
Symanzik equation to constrain the possible form of corrections to the classical potential 
in an arbitrary field background. 

For SU(2) Yang-Mills theory we found that the effective potential in a covariantly 
constant field background with T = \F^ v F^ va 7^ 0, Q = ^F^ u F^ va = is (suppressing 
the trace over the colour indices): 

Applying the Callan-Symanzik equation 



d d d 

^ + M dg-^ F dF. 



V eff =0 (2.78) 



we find that (3 = = — 3( -^ 2Af . These are properties of the Lagrangian, and do not 
depend on the particular background we evaluate it in; moreover to 1-loop order they 
are independent of the renormalization scheme. 

We now look for more general functions V that solve (12.781) . to see what possible 
corrections may appear in other field backgrounds. The equation (12.781) can be solved 
by a series of the form 

00 

V = J2^(9)F 2 log(gF/fi 2 y (2.79) 
where the a,(g) satisfy a set of coupled differential relations of the form 



ig ^i _ 2l a t + (i + l)a i+1 = (2.80) 
dg 

where we have used the relation /3 = 75 that we found above. 
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If we assume that to 1-loop order, the correction series in a particular background 
terminates at some order k, then we can integrate the relations (12.801 ) and impose that the 
function V reduces to the classical potential V = \F 2 plus corrections that are higher 
powers of g. We find 



a k = 

a k-i = Cig 2 

at _ 2 = c^-^i 

a 

a k - 3 = ••• (2.81) 

where we define the 1-loop (3 function (3(g) = ag 3 , a = — 3 ^n N - Thus, consistency 
with tree level fixes k = 2 and the value of C\, and subject to the assumptions above, the 
general effective potential for a (not necessarily constant) background with T ^ 0, Q = 
Ois 



V = ^ 2 + C 2 /F 2 -^F 2 log(^F 2 //i 4 ) 

The unfixed constant C 2 reflects the ability to shift the arbitrary renormalization scale 
fi, as well as the possible instability of the field background if C 2 is complex. Similar 
arguments constrain the form of V in an arbitrary background with Q ^ 0, which gives a 
generalization of the Euler-Heisenberg Lagrangian [BMS77 ]. Note in particular that the 
sign of the 1-loop contribution - and therefore the existence of the unstable perturbative 
vacuum - depends on the negative sign of (3(g). 
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Note that this method does not rely on knowledge of the precise the form of F" 
in 4-dimensional space-time, or in the internal (colour) space. Non-constant field con- 
figurations may have complicated derivative terms in their effective Lagrangian, but 
for configurations that satisfy our assumptions, the trace anomaly constrains the non- 
derivative terms to reduce essentially to the form of the constant field result obtained 
above. However, as noted above this does not allow us to reliably estimate the vacuum 
expectation value (F 2 ), because the 1-loop approximation still breaks down before we 
reach the dynamical scale A characteristic of confinement 6 . 

In section l2~3l we will turn this argument around, and use 1-loop anomalies to com- 
pute the effective superpotential of J\f = 1 supersymmetric Yang-Mills theory directly. 
The 1-loop anomaly calculation is exact in supersymmetric theories, which allows us to 
find the exact effective superpotential without needing to perform an explicit path inte- 
gral calculation around the vacuum field configuration. Indeed, the precise nature of the 
J\f = 1 vacuum is unknown, although we can compute some of its properties exactly. 

2.3 J\f = 1 supersymmetric gauge theories 

In a supersymmetric theory, the Lagrangian may contain terms of the form 



where the integral is over half of superspace, and W is the superpotential of the theory. 
It has dimension 3 and is a function of the chiral superfields $ 4 and not of their antichiral 

6 A more reliable estimate of (F 2 ) for QCD was made by Shifman et. al. ISV Z79 1 using charmonium 
sum rules. 




(2.83) 
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hermitian conjugates $j. The supersymmetric vacua of the theory are determined by the 
"F-term" constraints 



dW 



= 



(2.84) 



modulo complexified gauge transformations. In terms of the superpotential, the ordinary 
bosonic potential of the theory is given by 



where 0j are the lowest components of the chiral superfields $j and D a = £V |0,| 2 t a , 
where t a are the generators of the gauge group. 

There are two key results that allow us to compute the effective superpotential 
exactly in many supersymmetric theories: in a Wilsonian approach where we integrate 
over loop momenta down to a momentum cutoff, the superpotential only receives one- 
loop and non-perturbative corrections; and it is a holomorphic function of the chiral 
superfields and coupling constants. The meaning of these statements is somewhat sub- 
tle, and bears further explaining. 

Until now, we have considered the effective potential defined by the non-derivative 
terms in the generating functional of 1 -particle irreducible (1PI) diagrams of the the- 
ory that is obtained by integrating over the fluctuating fields. We found that in four- 
dimensional gauge theories this object receives contributions to all loop orders in per- 
turbation theory, corresponding to Feynman diagrams in the background field with arbi- 
trarily many internal loops. This remains true in a supersymmetric theory. Moreover, 
higher loop corrections will generically not be holomorphic. 

The Wilsonian approach to the effective action is to integrate over all loop momenta 
down to some cutoff scale; the resulting functional depends on the lower-momentum 




(2.85) 
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modes but has no dependence on momenta higher than the cutoff. If we integrate all 
the way to zero momentum we would recover the 1PI generating functional. In super- 
symmetric gauge theories, Shifman and Vainshtein HSV86I showed that the 2-loop and 
higher contributions are infrared effects; they only enter the Wilsonian effective action 
as the cutoff is taken to zero, and in computing matrix elements of Wilsonian quantities 
(averaging them over the external fields). For finite cutoff, the terms appearing in the 
Wilsonian effective action arise only from tree-level and 1-loop contributions. 

It is important to note that the parameters (fields, coupling constants) that appear in 
the Wilsonian effective action are not the physical quantities that would be measured in 
an experiment; indeed, the latter receive corrections to all orders. It would appear that 
the Wilsonian approach is missing the effects of the higher-loop contributions; as we 
saw in non- super symmetric Yang-Mills theory the higher loop corrections are vital for 
understanding the vacuum structure, because they dominate at low energies. 

The resolution, emphasized by HSV911IDS94L is that the all-loop, non-holomorphic 
1PI effective superpotential may be brought into the 1-loop, holomorphic Wilsonian 
form by a suitable (non-holomorphic, field- and coupling- dependent) change of vari- 
able. In other words, the 1PI effective superpotential is resummed into the Wilsonian 
form by this change of variable. This means that in supersymmetric theories the higher 
order corrections to the effective superpotential arising from the trace anomaly must all 
be related to the form of the 1-loop term, written in different variables. For example, in 
J\f = 1 supersymmetric Yang-Mills theory this is intimately related to the existence of 
the exact NSVZ /3-function [NSVZ83|, which has the form of a geometric series. 

Therefore, for supersymmetric theories we can confidently use the 1-loop Wilsonian 
effective potential to study the theory beyond the range where 1-loop perturbation theory 
naively breaks down, because we know that written in terms of physical quantities the 
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1-loop calculation sums the contributions to all loop orders. If in addition the non- 
perturbative corrections to the effective superpotential are calculable (by holomorphy 
and symmetry constraints, this is often the case), then we can obtain the exact effective 
superpotential, and by extension, exact results about the vacuum of the theory. The price 
is that to rewrite these exact Wilsonian results in terms of physical quantities one must 
undo the complicated change of variables. 



The effective superpotential for J\f = 1 Yang-Mills was constructed in [VY82J, by 
writing an effective Lagrangian whose symmetry transformations reproduced the correct 
1-loop anomalies. This is essentially the approach we used in earlier sections. 
The Lagrangian for A/" = 1 Yang-Mills theory is: 



where we have suppressed the gauge-fixing, ghost and auxilliary terms. In superfield 
notation this can be written as 



2.3.1 J\f = 1 Yang-Mills 



1 



_ F * vF ^a + Qpa^va + ^p a6 A 6 + . . . 



(2.86) 



4g 




(2.87) 



where we define 



S 



1 



TrW* 



32tt 2 



T 




+ i6 



(2.88) 
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S is the "gaugino bilinear superfield", whose lowest component is TrA 2 . In particular, 
S and r are both complex. 

The expansion of the composite superfield S in terms of component fields includes 
a term Tr(F^) 2 quadratic in the Yang-Mills field-strength tensors, which one might be 
tempted to identify with a scalar "glueball" operator of the Yang-Mills theory. How- 
ever, S cannot be interpreted as a dynamical glueball superfield, because the Yang-Mills 
field- strengths appear as auxilliary fields in S and are therefore non-dynamical HSS03H . 
The approach of studying the vacuum of J\f = 1 Yang-Mills theory by introducing a 
non-dynamical composite field is essentially the same approach we took in probing the 
Gross-Neveu model for the existence of a symmetry-breaking fermion condensate; here 
we are probing for a gaugino condensate, to which we associate the composite field S 
that includes the gaugino bilinear. In this sense, the effective superpotential W(S) we 
will obtain is part of a "minimal Lagrangian" that describes the symmetries and anoma- 
lies of the theory, but is not an effective Lagrangian for physical degrees of freedom. In 
particular, upon extremizing the effective superpotential W(S) we will obtain the value 
of the gaugino condensate in the vacua of J\f = 1 Yang-Mills. 

As before, the Callan-Symanzik equation constrains the form of corrections arising 
from the anomalous breaking of scale-invariance 7 : 



As we have seen in previous examples, it can be solved by a function of the form 



7 In Af — 1 Yang-Mills theory the trace anomaly is part of an anomaly multiplet that also includes 
the axial anomaly, and a superconformal anomaly. By supersymmetry, the constraints from the other 
anomalies are equivalent to that of the trace anomaly. 




(2.89) 




(2.90) 
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and we find 7 = 0, d = 8tt 2 , C 3 = = N, where (3{g) = -g^ to 1 loop. 

Therefore 

Wgff{S) = tS + C 2 S + NS 'log(4) 

= C 2 .S + A^log(S/A 3 ) (2.91) 

where we introduced the dynamical scale A via the running coupling relation 

r(/i) - log /i = log A (2.92) 

As in other examples, the constant C2 is not fixed by symmetries and may depend on 
the renormalization scheme. A value can be fixed following the approach of [CDSW02 1. 
Using an instanton calculation INSVZ 83 1, the value of the gaugino condensate can be 
obtained directly, giving rise to the value of the superpotential in the vacuum: 

W eff (A) = N(A 3N ) 1/N (2.93) 
The field S can be introduced by performing a Legendre transformation 




(2.94) 



Integrating out S recovers the previous expression (12.931) . If instead we integrate out C, 
then we recover the Veneziano- Yankielowicz superpotential 

W^S,A)=NS0og(^)-l) (2.95) 
which fixes the constant C2 = —N. 
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Since the field S here is complex, the F-term constraint = gives N distinct 
vacua (related by a phase, i.e. vacuum angle 9) 

(S) = e 2mk ^A 3 k = 0, . . . , N c - 1 (2.96) 

Furthermore, as noted in the previous section, this Wilsonian effective superpotential 
does not receive corrections beyond one loop. Therefore the vacuum expectation value 
(S) oc (TrW a W a ) is exact, and the N vacua of J\f = 1 supersymmetric SU(N) Yang- 
Mills theory have a non- vanishing gaugino condensate. 

Note that the Callan-Symanzik anomaly calculation does not assume a particular 
form of the background gauge field configuration. A covariantly constant background 
field strength was considered in [Kay83|, generalizing the Yang-Mills calculations 
reviewed in section 12.2.21 As in the non- supersymmetric case, a constant background 
field strength causes the vacuum energy to decrease, but there is still an instability at the 
1-loop level 8 . A field theoretical derivation of the Veneziano-Yankielowicz superpoten- 
tial is not known - this would amount to knowing the field configuration in the J\f = 1 
Yang-Mills vacuum and integrating over the fluctuations around this background. 

2.4 J\f = 1 theories with matter 

One of the starting-points for the recent work on J\f = 1 gauge theories with adjoint 
matter was the conjecture [DV02a, DV02"bl lbV02cl that the effective superpotential is 
computed by an associated bosonic large- N matrix integral, which may be evaluated by 
counting planar diagrams. This conjecture comes from string theory, and follows a chain 
of reasoning that is the culmination of extensive research on the relationship between 
string theory and gauge theories. 

8 This is not surprising since this field configuration is not supersymmetric. 

42 



The steps in the conjecture can be summarized as follows: type II string theory on 
certain Calabi-Yau manifolds ("generalized conifolds") is known to reduce to Af = 1 
Yang-Mills theories in a limit that decouples gravity; at low energies these geometrical 
spaces undergo a "geometric transition", where a cycle in the geometry shrinks to zero 
size and is replaced by a different cycle of finite size. This is a geometrical analogue 
to confinement of the Yang-Mills theory at low energies. If we instead consider B-type 
topological strings on these spaces, the topological string amplitudes reproduce the F- 
terms (superpotential) of the corresponding gauge theory. Therefore, after the geometric 
transition they should give us the gauge theory effective superpotential. However, the 
path integral of the topological B-model on these spaces reduces to a large N matrix 
integral. Following the chain of arguments, the effective superpotential of Af = 1 Yang- 
Mills theories should reduce to a large N matrix integral. Thus, string theory provided 
an entirely unexpected computational tool for studying the effective superpotential of 
Af = 1 gauge theories with matter. 

In practical terms, we can illustrate the technique as follows. Suppose we start with 
a SU(N C ) gauge theory with Af = 1 supersymmetry and a chiral superfield $ in the 
adjoint representation, with a tree-level superpotential that contains a mass term and 
cubic self-interaction: 



in terms of the gaugino bilinear S, receives contributions from two sources: Veneziano- 
Yankielowicz terms arising from the strongly-coupled dynamics of the gauge field, and 
contributions from the matter field $. According to the conjecture, the only contri- 
butions of the matter field $ to the effective superpotential come from the planar $ 
diagrams of the theory (even at finite iV c ) where we insert the external S field once into 




(2.97) 



String theory suggests that the effective superpotential of this theory W e jj{S), written 
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each of the index loops of the $ diagrams. Furthermore the effective superpotential has 
no dependence on the internal loop momenta of the diagrams ! 

The meaning of this result is that the superpotential for such theories is an essentially 
combinatorial object, depending only on the counting of ribbon diagrams with planar 
topology. It has been known for a long time that these planar diagrams are counted by 
a zero-dimensional matrix integral |BIPZ78 1, and we can often evaluate the free energy 
of this "matrix model" exactly. 

We saw in the previous section that in non-supersymmetric field theories the need 
to integrate over loop momenta was a serious complication for extending the compu- 
tation of the effective action to higher orders. What is the field theory process that 
removes the contribution of loop integrals when supersymmetry is present? As in non- 
supersymmetric theories, we can understand the field theory results in two ways: using 
anomalies [CDSW02| and by evaluating the path integral llDGL + 03l. We will summa- 
rize the results of these papers, and refer to the original papers for the details. 

The technique of using anomalous symmetries to solve for the effective superpoten- 
tial has been extended to a large class of J\f = 1 theories HCDSW021 ISei03l IBIN+03L 
where the relevant anomalies are of generalized Konishi type. This approach relies on 
the fact that the set of chiral primary fields - those that can enter the effective superpo- 
tential - are closed under addition and operator product, up to terms that vanish when 
evaluated in a supersymmetric vacuum; in other words the chiral primary fields generate 
a ring structure, the chiral ring. Moreover, elements of the chiral ring are independent 
of position, so the chiral ring is a global structure. 

Using the properties of the chiral ring, it was shown that the (anomalous) symme- 
tries of the theory (particularly the generalized Konishi anomalies) restrict the possible 
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superpotential contributions to the planar diagrams with insertions of S 9 . Then, the 
Ward identities associated to the generalized Konishi anomalies are shown to be equiv- 
alent to the loop equations of the matrix model, which are Dyson-Schwinger equations 
for the correlation functions, and which can be solved using matrix integral techniques 
to determine the effective superpotential exactly. 

A complementary field theory approach llDGL + 03l used the background field 
method to study J\f = 1 gauge theories. They showed that as a consequence of sym- 
metries, it is again only the planar diagrams of the gauge theory that can contribute to 
the effective superpotential, and moreover supersymmetry implies that after the loop 
diagrams are summed in the Schwinger formalism, the loop momentum dependence in 
the diagram sum exactly cancels between bosonic and fermionic contributions. Since 
there is no remaining dependence on loop momenta, the resulting effective superpoten- 
tial reduces to the zero-dimensional matrix model calculation. A key feature seen in this 
approach is that the individual gauge theory loop diagrams do depend on loop momenta, 
but after summing over all diagrams the momentum dependence exactly cancels. 

There are several remarkable consequences of these results. In many cases the asso- 
ciated matrix integral can be directly solved (corresponding to summing the Feynman 
diagram expansion to all orders). However, in more complicated examples where the 
diagram series cannot easily be summed using known techniques, a perturbative expan- 
sion of the ribbon diagrams (up to some order in the number of index loops) gives a 
perturbative expansion of the effective superpotential W(S), which upon extremization 
generates an expansion of the vacuum gluino condensate (S) ~ (AA) as a sum of frac- 
tional instanton contributions. As emphasized in IDV02cl . and as we have seen in other 
examples above, the perturbative loop expansion of the gauge theory in terms of an 

9 By contrast to the trace anomaly, the generalized Konishi anomalies contribute to all orders of per- 
turbation theory, although in a simple, and often summable way. 
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appropriate choice of composite operator yields non-perturbative information about the 
vacuum. 

These results have been checked and extended in a large number of papers, and the 
deeper consequences for the quantum structure of gauge theories are still being explored. 
In the remainder of this thesis we will discuss our contributions to this area of research. 
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Chapter 3 

Effective Superpotentials from 
Geometry 

In this chapter we first review how gauge theories with J\f = 1 supersymmetry may 
be obtained from string theory, and how string theory provides new tools for analyzing 
their low energy structure. The simplest examples have U(N) or SU(N) gauge group 
with matter in the adjoint and fundamental representation. In llACH + 03l we extended 
the analysis to SO(N) and Sp(N) gauge groups with adjoint matter, and in |KW03| we 
showed that a careful consideration of UV divergences requires the inclusion of a maxi- 
mal number of fundamental matter fields in order to regulate those divergences. We then 
studied in detail the structure of U(N) and SU(N) theories with adjoint and fundamental 
matter and developed simple, general formulae for the effective superpotentials, which 
reduce in special cases to previously known results. Other J\f = 1 theories have been 
treated in the literature including theories that do not arise from soft supersymmetry 
breaking of an J\f = 2 theory HLLT04I . 

3.1 Geometric engineering of gauge theories 

We begin by reviewing the construction from string theory of a softly broken J\f = 2 
gauge theory with SUISOISp gauge group llCv^OWTllEOT()l[lDV02a1l . Consider 
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type IIB string theory compactified on the non-compact A\ fibration 

u 2 + v 2 + w 2 + W'(x) 2 = 0, (3.1) 

where W(x) is a degree n + 1 polynomial, which will later be related to the tree level 
superpotential for the adjoint chiral superfield $. This fibration has singularities at the 
critical points of W(x). In the neighborhood of those singularities, we can introduce 
the coordinate x' = W'(x). Then it is easy to see that the singularities are all conifold 
singularities. 

This generalized conifold can be de-singularized in two ways: it can be resolved or 
it can be deformed. The resolution is given by the surface 

u + iv w + iW'(x) \ I Ai 
— w + iW'(x) u — iv J \ A 2 

in C 4 x P 1 . In this geometry each singular point is replaced by a P 1 . These P^s are 
disjoint, holomorphic, have the same volume and are homologically equivalent. The 
latter property can be seen by making use of the fibration structure away from W'(x) = 
0. This A 1 fibration over the x plane induces a fibration of non-holomorphic S 2, s over 
the x plane. This S 2 cannot shrink to zero size as one approaches a critical point of W 
in the x plane, but it becomes the holomorphic P 1 of the resolution. 

We can now construct a softly broken J\f = 2 U(N) gauge theory with tree level 
superpotential W(x) by wrapping N D5-branes around the S 2 . The adjoint chiral 
superfield $ parameterizes the normal deformations of the D-branes, and since these 
deformations are obstructed in the Calabi-Yau geometry there is a superpotential for $, 
which is identified with the function W(x) that describes the nontrivial A\ fibration of 
the generalized conifold IBDLROOl lKKLMOOII . This is an UV definition of the theory; 
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in fact it describes quantum gravity coupled to the gauge theory, because the excita- 
tion spectrum also includes the closed strings that propagate away from the D-branes, 
and which give rise to gravitons in the particle spectrum. The bulk modes and massive 
open string modes can be decoupled by taking the 't Hooft limit N ^ oo, g s ^ 0, 
A = g s N = const., which leaves only the lowest open string modes, the gauge and 
matter fields. 

A classical supersymmetric vacuum of the gauge theory is obtained by minimizing 
the volume of the D5-branes. This amounts to distributing a collection of Ni D5-branes 
over the n minimal- volume holomorphic P x 's at the critical points of W. The U(N) 
gauge symmetry is then spontaneously broken to U(Nx) x •■ ■ x Z7(iV Tl _ 1 ). SU(N) 
gauge group can be treated by decoupling the overall U(l) C U(N) trace, which is a 
free theory. 

If we flow this ultraviolet theory to the infrared (low energies), there will be a 
confinement transition. In string theory this is described by a "geometric transition" 
in which the resolved conifold geometry with wrapped D5-branes is replaced by a 
deformed conifold geometry HVafOlH 

u 2 + v 2 + w 2 + W'(x) 2 - f(x) = 0, (3.3) 

where f(x) is a polynomial of degree n — 1. For a reasonably small f(x), each critical 
point of W'(x) is replaced by two simple zeros of W'(x) 2 — f(x). This means that 
each Fj is replaced by a 3-sphere A, with 3-form RR-flux H through it, equal to the 
amount of D5-brane charge on the Fj. After the geometric transition there are no more 
D-branes, so there are only closed strings in the spectrum. 
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The coefficients in f(x) are normalizable modes that are localized close to the tip of 
the conifold. The coefficients in fix) are determined by the periods 



These periods Si are to be identified with the gaugino bilinear superfields of the gauge 
theory. There are non-compact 3-cycles Bi that are dual to the A^ The periods of the 
B-cycles are 



where JF is the prepotential of the Calabi-Yau geometry. One needs to introduce a cutoff 
in order to make these periods finite; we will discuss the physical meaning of this cutoff 
in section l3~L4l 

The flux through the cycles Ai is determined in terms of the RR-charges of the D- 
brane configuration 




(3.4) 




(3.5) 




(3.6) 



and the flux through the cycles Bi is given in terms of the coupling constants 




(3.7) 



The effective superpotential W e jj{Si) is then given by the flux superpotential HTVOOl 



BB96l lGVW00l fPS96l 




(3.8) 
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Figure 3.1: The complex curve that results from projecting the Calabi-Yau to the 
base of the S 2 fibration. It is a branched double cover of the complex plane, where 
the cuts are the projections of the S 3 cycles of the Calabi-Yau. The A contours are 
compact cycles, and the B contours = BT + Bf are non-compact and run from 
a point at infinity on the lower sheet, through the i th cut to the point at infinity on 
the upper sheet. For later convenience the B contours have been regularized by a 
cutoff A . 



Using the expressions for the periods and the fluxes, we get 



In evaluating these period integrals, the u and v integrals can be performed trivially 
(the A and B cycles have the form of an S 2 fibration over lines in the complex plane, 
see figure l3~TT) . and the period integrals of the complex 3-dimensional Calabi-Yau (13.31) 
can be reduced to the period integrals of a complex curve 




(3.9) 



y 2 = W{x) 2 - f(x) 



(3.10) 
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The holomorphic 3-form fi, the periods of which define the effective superpotential, 
reduces to the meromorphic 1-form y dx on the curve. The function f(x), and therefore 
the curve itself, is fixed by a requirement of extremality, in a sense that will be made 
precise. This curve is central to the construction of the gauge theory effective superpo- 
tentials, and we will rederive and study it from several points of view in the following 
chapters. 

In order to study SO or Sp gauge theory, we can consider an orientifold of the 
previous geometry 1 . Since we started with a type IIB theory on a Calabi-Yau, we have 
to combine the worldsheet orientation reversal with a holomorphic involution of the 
Calabi-Yau (an anti-holomorphic involution would be appropriate for the IIA theory). 
Furthermore we want to fix one of the P x 's and act freely on the rest of the Calabi Yau 
geometry. This can be done if W(x) is an even polynomial of order 2n. In terms of the 
fibration structure of the Calabi-Yau, this means that the critical points of W (x) come 
in pairs (— Xi, Xi) and one critical point is fixed at x = 0. Then 

(u, v, w, x, Ai, A 2 ) i-> (— u, —v, —w, -x, Ai, A 2 ) (3.11) 

is a holomorphic involution of the geometry (I3.2I) . which leaves only the P 1 at u = v = 
w = x = fixed. In the string theory this means that there is an 05-plane wrapping this 
P 1 in the Calabi-Yau geometry. 

There are essentially two choices of 05-plane with which we can wrap the fixed P 1 . 
They are distinguished by a different choice of worldsheet action and carry RR 5-form 
charge of ±1 (the RR charge of an -plane is ±2 P_5 in conventions where we count 
the charge of N/2 D-branes but not their N/2 images). The orientifold contribution 
to the RR charge of objects wrapping the P 1 will cause a shift in the coefficient No in 

'Orientifolds were discussed in the A-model in IS V00IIEOT0TllAAHV02IIFO03l . while the discussion 
of |Gom02| is more closely related to the B-model which is our interest here. 
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the flux-generated superpotential on the deformed Calabi-Yau geometry, as explained 
below. 

Now we can construct a softly broken J\f = 2 SO(N) or Sp(N/2) gauge theory 
with tree level superpotential W(x) by wrapping N D5-branes around the S 2 and then 
performing the orientifold. The gauge symmetry is again broken SO(N) i— > SO(N ) x 
U{Nx) x • ■ ■ x C/(iV n _i) or Sp{N/2) i-> Sp(N /2) x C/(JVi) x • • • x C/(iV n _i) respectively, 
with iV = iV + 2A^x + h 2iV n _i. 

At low energies the geometric transition again produces the deformed conifold 
geometry [VafOl | 



where f(x) is now an even polynomial of degree 2n — 2. Such a polynomial represents 
the most general normalizable deformation of the singular conifold that still respects the 
holomorphic involution (13 .lib . The orientifold acts on one 3-sphere A as the antipodal 
map, while the other 3-spheres are mapped to each other in pairs Ai and A-*. Note that 
there is no orientifold fixed plane anymore. 

The 3-form RR-flux H through each 3-sphere is equal to the amount of D5-brane 
and 05-plane charge on the ¥} before the transition. 



and the flux through the cycles B{ is again given in terms of the coupling constants 



u 2 + v 2 + w 2 + W'{x) 2 - f{x) = 0, 



(3.12) 




(3.13) 




(3.14) 
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Since there is no orientifold fixed plane, there are no contributions to the effective super- 
potential for the gaugino condensate from unoriented closed strings IAA HV02I . In the 
flux superpotential 



the integral is now taken only over half of the covering space of the orientifold. Using 
the expressions for the periods and the fluxes and taking into account the orientifold 
projection, we get 



This result could also have been computed on the open string side before the tran- 
sition. On the open string side there is no flux through any 3-cycles, so there is no 
contribution to the superpotential due to closed oriented strings. But there are two kinds 
of other contributions to the effective superpotential: the open string contributions (disk 
diagrams) and the contributions due to closed unoriented strings at the orientifold fixed 
plane (MP 2 diagrams). The contribution due to the open strings is the equal to one half 
that of the theory without the orientifold, i.e., it is 




(3.15) 




(3.16) 




The contribution due to the unoriented closed strings then must be 



WtfiSi) = W ej] {Si) - W% f {Si) = ± 



dSo' 



(3.18) 



We will confirm this result in a matrix model computation in chapter[5] 
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3.1.1 Computing the superpotential 

Consider pure J\f = 2 Yang-Mills theory broken to J\f = 1 via a tree-level superpotential 
of the form: 



which is single- valued on the genus n — 1 Riemann surface y 2 = W'(x) 2 + f n -i(x) 
(the "A/" = 1 curve") that we encountered in the previous section. In section 13.21 we 
will rederive this curve by factorizing the Seiberg-Witten curve of the associated J\f = 2 
theory obtained when Wf ree = 0, and discarding the repeated roots of the curve that 
correspond to condensed monopoles. 

The compact A -periods of the curve yield the gaugino bilinear superfields, Si, while 
the non-compact B -periods, IT yield the derivatives of the free energy Choose the 
branches of the square root so that on the first sheet u(x) vanishes in the classical limit 
fn-i 0; therefore on the second sheet u(x) — ► W'(x). 

In this chapter we will focus on the maximally-confining phase of the theory (the 
vacua with classically unbroken gauge group U(N)), for which the resolvent degener- 
ates: 



y{x) = V(W'(x)) 2 + fn-^dx = G n _i(x) y/(x-c) 2 - fi 2 dx, (3.21) 




(3.19) 



The effective superpotential W e -a(S) may be computed in terms of periods of the differ- 
ential form ("resolvent"): 
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for some polynomial, G n -i(x) of degree (n — 1). For U(N) theories, it is convenient to 
use the freedom to shift x so as to set c = 0; this is not allowed for SU(N), for which the 
center of the cut is not a free parameter, but the SU(N) results may be obtained from the 
U(N) at the end of the calculation by decoupling the overall U(l) trace (we will come 
back to this point later). The gaugino bilinear is then given by: 

S = J_ / u ( x ) = ±_L / y ( x ) = ± J_ T G^x) ^x 2 - fi 2 dx (3.22) 
2tt2 J a Am J A 2m 

where the sign depends on the orientation of the contour. The £>-period is given by 
integrating along a contour from infinity on the second sheet, through the cut to infinity 
on the first sheet, see figure l3~Tl The logarithmic divergence of this integral needs to be 
regularized, and this is usually done by a introducing a UV cut-off: 

x + =A i-Ao 



Ub — u — lu = — G n -i(x) y x 2 — \i 2 dx , (3.23) 

J B Jx-=Ao J /j, 

where x_ and x + denote the values of x on the lower and upper sheets respectively. The 
effective superpotential is then given by: 

W eS = NUb + NW(A ) + tS (3.24) 

where r is the bare gauge coupling, and the second term is added to cancel the con- 
tribution from the upper limit of the integral in Ub- As we saw in section 12.3.11 the 
effect of the r term is to combine with the log-divergent piece of Ub to give the (finite) 
dynamical scale A of the theory. 

In computing the effective superpotential by this method, the approach initially taken 
in the literature was to send A — > oo, causing its effects to decouple from the theory. 
However, we will obtain physical insight into the nature of the computation by keeping 
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the cut-off finite. We will henceforth take the cut-off A to be large but finite, and 
investigate the effects on the low-energy gauge theory; this amounts to keeping the 
O(l/A ) terms in Tl B and subsequent calculations. 

3.1.2 Example: U{2) 

Before analyzing the general case, consider the simplest example of U(2) with a tree- 
level mass: W = |mTr<3> 2 . The effective 1-form is 

y(x) = m\j x 2 — jj? (3.25) 

which is single-valued on a two-sheeted Riemann surface with a cut between x = ±fi. 
The gaugino bilinear is given by the A-period: 




(3.26) 



(3.27) 

where the integral is evaluated using hyperbolic functions and the two branches come 
from sinh(:r) = ±A/cosh 2 (:r) — 1 (this amounts to a choice of contour, i.e. integrating to 
the point above A on one of the two sheets). As mentioned in the previous section, the 
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role of r in (13.241) is to replace the N log(mAQ) term in lie by the finite scale iV log(A 3 
This may be implemented in practice by setting r = N log(^?) in (13 .241 ) . 
We find 



„. . (( , , , ,Sa S 2 2S 3 5S 4 US 5 

W = N [ S 1 - lo 



A 3 V mkl 2(mA 2 ) 2 3{mA 2 ) 3 4(mA 2 ,) 4 " 

(3.28) 

Therefore in the limit A — > oo (equivalently, keeping A finite and considering energies 
m « A ) the infinite correction series tends to zero and the effective superpotential 
(13.281) reduces to the usual Veneziano-Yankielowicz superpotential. 

The form of the series (13.281) is the same as that obtained for U(2), Nf = 4, with A 
identified with the quark mass. The known formula for W(S) with tree-level superpo- 
tential W = |mTr$ 2 + Y^fJi ^QiQ 1 + QMQ { is HACFH()3bllBrN+03ll 



W(S) = N C S(1 - log(^f)) - N f Sh ^ 



-N f S g + - log(i±4EM)) (3 , 9) 

with a = l/(m/i 2 ) (we will derive this expression in section I3.1.4I) . Setting A^ c = 
2, Nj = 4, ii = A and performing the series expansion, we recover the expression in 
(13.281) . We will show in Section 13.1.41 that this feature remains true for general A^ c and 
W(&), and the corrections obtained by keeping the cut-off dependence in the period 
integral indeed have the physical interpretation of Nf = 2N C massive quark superfields, 
which serve to regularize the divergences of the calculation. 

Choosing the other branch of I1 B we obtain the negative of (I3.28I) . This branch 
describes a Higgs branch IBI N + 03j . where the gauge symmetry is broken by giving a 
vev to the scalar component of the quark superfields (an arbitrary Higgs vacuum can 
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be obtained by writing W = tS + Y^i=i and choosing the branch of lis termwise, 
i.e. for each period integral we choose whether to integrate along a contour on the first 
or second sheet). 

If instead of U(N) gauge theory we considered SU(N), the foregoing discussion 
would be modified by the need to ensure "quantum tracelessness" of the vacuum, i.e. that 
(tti) — 0. This may be achieved by taking the tree-level superpotential W = |mTr$ 2 + 
ATr$ and proceeding with the above analysis, treating A as a Lagrange multiplier to 
enforce (m) = (Tr$) = 0. Instead of repeating the calculation for SU(2), we will defer 
until later when we consider the general U(N) and SU(N) cases. 

3.1.3 Evaluation of the period integral for general W 

The period integrals, (13.221) and (13.231) . are elementary but one can obtain a simple closed 
form in terms of Wf ree . This can be evaluated and gives a combinatorial formula for the 
moduli Uk which can be compared to other techniques. Make the change of variables 2 : 

x = + r 1 ), (3-30) 

and define series expansions: 

1 n+l 

W(-n(t + r 1 )) = b + Yl b *(Z k + (3.31) 

k=l 

1 n 

+ r 1 )) = c + Yl c ^(^ + (3.32) 

k=l 

2 We again assume that x has been centered on the cut. 
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Note that the series take this form because of the symmetry of (13.301) under £ — > £ 1 . 
Under this change of variables the integrand may be written: 

i|u(f - r')G„-i(if«({ + r 1 )) = C-iWVi' - /i 2 , (3.33) 

= ^f+Mfr* <334) 

= W'(a;) + C^- 1 ) . (3.35) 

The left-hand side is manifestly odd under £ — > l/£, while the right-hand side shows 
that all the non-negative powers in the £-expansion are given by (13.321) . It therefore 
follows that under the change of variables, one has 

n 

s/{W'{x)Y + / n _!(x) = G n ^{x) y/{tf - ^) = c * (£* ~ r*) • (3-36) 



k=l 



Note in particular that the left-hand side of (13.331) is manifestly odd under £ — > l/£, 
therefore c = in (13.32b . 

Define [...]_ to mean: discard all the non-negative powers of £ in [...]. We may 
then write the last equation as: 



. (3.37) 



One can now easily perform the integrals (13.221) and (13.23b . The former is simply given 
by taking £ = e %e for < 9 < it, and it picks out the £-residue: 



S = | ci (3.38) 
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To perform the second integral first note that: 



d_ 



i 



2^(i-r 2 ) 



(3.39) 



and therefore: 



^(W'(X)) 2 + f n -l{x)dx 



-fici log(0 + W(z) -2 
-/iCi log(f) + b 

n 



(3.40) 



(3.41) 



k=l 



where (£ + £ To obtain II, we must evaluate this between £ = 1 and £ = £ > 

where 

An / / / a \ 2\ 

' » ' (3.42) 



This yields: 



n = bo + jic logKo) - (»'(A )-2[»'( 



(3.43) 



where the definite integral has been evaluated using (13.401) at £ = A and using (13.411) at 

£ = i. 

In the limit of large A the last term in (13.431 ) vanishes since it only involves negative 
powers of £ ~ A Taking this limit, and using ( 13.381) one obtains: 



II = b + 2Slog(^) - W(A ). 



(3.44) 
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Therefore 



Nb + 2NS log (— 
V fj, 



) 



(3.45) 



We will show in section IT. 1.51 that for general Wt ree (&), (13.451) can be extremized with 
respect to S by taking \i = 2 A, and we find the previously known result ICIV01I 



where we have evaluated the coefficients 6 in the series expansion (13.311) . 

In the previous section we discussed the geometric engineering of this gauge theory 
from string theory, which involved D-branes wrapped on cycles of a Calabi-Yau. From 
the string theory perspective it is tempting to also interpret the cut-off of the period 
contour in terms of branes. That is, it is really only physically natural to terminate the 
period integral on another brane. Since D-branes carry gauge fields, having a stack of M 
branes at the point A would mean that one started with a larger (product) gauge group 
and that the original SU(N) theory is actually coupled to M bi-fundamental matter mul- 
tiplets with a (gauged) SU(M) "flavor" group (see [Hof03 1 for an analysis of this gauge 
theory). However, when the second set of branes become non-compact, their associated 
gauge coupling tends to zero, and the SU(M) gauge factor becomes a global SU(M) 
flavor symmetry. Thus, string theory suggests that keeping the UV cut-off terms should 
yield the superpotential associated with the coupling to fundamental matter multiplets. 
This is indeed what we find in explicit calculations. 

If one also recalls that the canonical form of the £>-period integral, (I3.23I) . involves 
an integral from the lower to the upper sheet of the Riemann surface, then this extra 




(3.46) 
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term may be thought of arising from iV c branes (or anti-branes) at each limit. Thus 
one can also extract the results for Nf = N c by regulating the upper and lower limits 
independently. We will develop and extend this observation in the next section. 



3.1.4 UV cut-off as regularization by Nf = 2N C fundamental 



As mentioned in section U.l.ll the effective superpotential for the Nf = theory (in a 
maximally confining vacuum) is given by IICIVO 1 1 lr~W03bl ICSW031 



where the integral is formally divergent and is usually cut off at a point A . We will ver- 
ify in section l5~2l that introducing Nf fundamentals gives the (again formally divergent) 
contribution IIDV()2cllACF~H03bll 



However, when Nf = 2N C , the contours combine and the integration domains are 
now finite, so the divergence of the integrals have been regularized. When all m; are 
equal we may write = A and we can explicitly see the role of the 2 N c fundamental 
fields in implementing the cut-off of the Nf = integral: they act as regulators for the 
UV divergences of the calculation, by removing the short-distance divergences of the 
calculation. Physically, the gauge theory with an adjoint chiral superfield and Nf = 
2N C fundamentals has vanishing beta function in the limit when all of the fields are 
effectively massless, i.e. at energy scales much greater than their mass. Thus, the theory 
has a nontrivial UV conformal fixed point, and is free from short-distance singularities. 



quarks 




(3.47) 




(3.48) 
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In terms of the additional microscopic degrees of freedom we are forced to add, the 
tree level superpotential of the gauge theory is modified: 

W t ree{®) - W tree {<S>) + ^ A W + Q W (3.49) 

i=i 

where Q l are the new "quark" superfields, and Q % are their conjugate antiquaries, and 
we have normalized the coefficient of the Yukawa interaction to 1 (the Yukawa coupling 
can be absorbed into the mass parameters m; by redefining the fields, since we are not 
interested in the kinetic terms). In section I5T21 we show how the combinatorics of the 
Feynman diagrams involving the new quark fields combine to subtract the short-distance 
divergences of the theory without quarks. 

As we have seen in the example of U(2), when A is taken to be large but finite, 
it gives finite (but small) corrections to the expression for the effective superpotential 
W(S). Therefore, the vacuum expectation value for the gaugino bilinears (Si) will be 
perturbed from that of the theory we started with (J\f = 1 Yang-Mills theory with a 
massive adjoint and no fundamental matter). In other words, in terms of the J\f = 1 
curve (13.101) . the presence of the cut-off at a finite distance from the cuts cause the size 
and center of the cuts to be perturbed. Because of this deformation, it will turn out that 
this J\f = 1 curve cannot be obtained by factorizing the SW curve of pure J\f = 2 Yang 
Mills. 

Therefore, in regularizing the Nf = theory by imposing a finite cut-off on the 
divergent integral, we have gone off-shell (i.e. the vacua of this theory do not solve the 
equations of motion of the Nf = theory). Physically, this is because the presence of 
the cutoff is equivalent to introducing new physical degrees of freedom that contribute 
to the gaugino condensates. This amounts to embedding the Nf = theory in a larger 
theory with Nf = 2N C massive quark flavors; it is only in the limit of infinite quark 
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mass (infinite cut-off) that the effects of the quarks on the vacuum structure of the theory 
decouple and we approach the on- shell vacua of the Nf = theory. 

In practice we can think of the effective superpotential for < Nf < 2N C fun- 
damentals (as computed using the technique of [CIV01 ] discussed here, also using the 
matrix model discussed in chapter|5]), as always being generated by the UV-finite theory 
with 2N C fundamental fields, with masses that are either kept finite or which are taken 
to infinity at the end of the calculation and decouple from the theory. In other words, if 
we have Nf fundamental fields of finite mass, then the remaining 2N C — Nf are of mass 
A ^> m. Therefore: 



poo J poo poo 

W eff ~ -2N C / uj + J^ uj + (2N c -N f ) 

" J H ,- = i J mi J A 



N f , 

oo poo \ ' / poo poo 



-2nA / u- / u 



UJ — UJ 
H J Ao / ~\ V^Ao Jmi 

fA N f pA 

2N C u + J^ v ( 3 - 5 °) 

J u „■_ i Jmi 



i=l " 



and all integrals are finite. 

We can then decouple the quarks of mass A by taking A — > oo, and using the 
results of section l3T~3l we find the following expression for W e ff : 
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W eff = N c (b + fia M^)) + 

Nf 

= ^ + ^, og ( (^)^na^) j + g| H/K(m;))1 _ +TS 

+ X^[H/(eW)]_ + r5 

/ r ) 2N c -N f '\2N c N f / I 7— 

+ (3-51) 



1=1 



where we used the scale-matching relation A 2Nc = A 2Nc ~ N f Yli m i- Using the defini- 
tions (13.311) and writing explicit expressions for the coefficients b k , this can be written 
as: 



2i » 92i (^ 2i 



/ o^Nc-Nf \2N C N f / : " \ 

N f n+1 , k / l \ 

+EEf(f) E W" a <« 2 > 

8=1 fc=l Z=|fc/2J+1 \ t / 
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An explicit general expression for S = S(gk, n) can be similarly obtained, but we will 
not need it here. 



3.1.5 Extremizing the superpotential 

In order to find the physical vacua, we need to extremize (13.521) with respect to S. This 
will fix /Li, the size of the cut, and give the vacuum superpotential in terms of physical 
quantities. Varying with respect to S, this can be achieved by setting 

l = ° 



, n2N c ~N f \2N c N f / I — , 



F i-i 




i.e. 



N f 1 + Jl-(^) 2 \ 

f Nc = {2K) 2N C JJ V _ (3 55) 

Thus, the logarithmic term in (13.521) does not contribute in the vacuum, and the extremal 



superpotential is found by solving (13.551) to find (//) . Note that when Nf — the solution 
to (13.551) is given by taking \i — 2A = 2A, as claimed in section l3~T31 

When all quark masses are taken equal, rrii = m, (13.551 ) can be written in the sim- 
plified form 

, K2N B /N f _ ^A^Y c/Nf + -^(4A 2 ) 2N ^ N f = (3.56) 

Am 2 

Note that this condition is polynomial in ji 2 when N c is a multiple of Nf. 
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3.1.6 Examples 



We turn now to some other examples (in all cases the quark masses are set equal for 
simplicity). 

Quadratic tree-level superpotential 

The simplest tree-level superpotential of the form (13.491) contains a mass term for the 
adjoint chiral superfield $: 



We consider arbitrary values of N c and Nf. The gaugino bilinear takes the simple form 




(3.57) 



S 




M 



(3.58) 



and we can eliminate \i from W e jj{m, M, A, //) to write the effective superpotential in 
terms of the physical parameters and gaugino bilinear: 



N c (S + Slog( 



M N C ~ K 2N C 



))+N f Slog( 



1 + VI - 45a 
2 



eff 



l-2Sa + VI - 45a 



N C S(1 + log( 




1 + VI - ASa 
2 



(3.59) 
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where a = j^z- This is the previously claimed result ( 13.291) . first obtained by 
[ACFH03b|. For the special case iV c = 2, Nf = 1 the extremization condition (13.561 ) 
becomes 



// - (4A V) 2 + 




(4A 2 ) 4 = 



(3.60) 



■A 



S 2 A 6 + SA 12 a = 



(3.61) 



where A 3 = MA 2 is the scale of the theory below the mass M of the adjoint. Excluding 
the unphysical solution S = (which would correspond to a vacuum with unbroken 
chiral symmetry, and can be ruled out on general grounds [CDSW02|), there are three 
remaining solutions. Taking the limit of infinite quark mass, a — >• 0, (13.611 ) degenerates 
further: 



i.e. two solutions S = are unphysical, and the two physical solutions are S = ±A 3 . 
At energies much lower than the mass M of the adjoint field $, the theory is described 
by M = 1 SU(2) Yang-Mills, and we indeed obtain the correct value of the gaugino 
condensates (12.961 ) of the Veneziano-Yankielowicz superpotential. 

Keeping the mass of the fundamental fields finite gives a series of corrections to the 
pure J\f = 1 result: 



S 2 (S 2 - A 6 



) = 



(3.62) 
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A 3 - i«A 6 - f a 2 A 9 - i« 3 A 12 - j|a 4 A 15 + . . . 
(S) = { 2 8 2 128 (3.63) 

-A 3 - i«A 6 + fa 2 A 9 - ±a 3 A 12 + ±fa 4 A 15 + ... 

2A 3 - iaA 6 - ia 2 A 9 - VA 12 - §a 4 A 15 + . . . 
Wi ow = { 2 4 4 64 (3.64) 

-2A 3 - ±aA 6 + ±a 2 A 9 - ±a 3 A 12 + |a 4 A 15 + . . . 

This result agrees with that of [ ACFH03b | (although they only explicitly considered one 
of the two vacua). It shows clearly how the presence of the finite-mass quarks perturbs 
the vacua of the theory away from their Nf = values. 

Equation (13.611 ) encodes the exact form of the effective superpotential of this theory. 
In this case a closed-form expression for (S) and Wi ow could also be obtained since the 
cubic branch of equation (13.611) may be solved explicitly; for higher-rank gauge groups 
the polynomial will be of degree 2N C — 1 in S, and can always at least be evaluated as a 
series expansion to any desired order. 

Arbitrary tree-level superpotential with Nf = N c 

In this example the extremization constraint (13.561) becomes quadratic in /i 2 , and can be 
trivially solved for arbitrary tree-level superpotential W / ($): 



// - (4Am - 4A 2 )/i 2 = 0, (3.65) 

so /i 2 = 4(Am — A 2 ) (the solution /i 2 = is again unphysical). 

When the cut in the J\f = 1 curve is centered away from the origin, centering the 
coordinate axes on the cut introduces a corresponding shift in the quark masses, m i— > 
m + c. In the following section we will see that factorizing the U(N C ) Seiberg-Witten 
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curve for Nf = N c fixes c = A, hence /i 2 = 4Am, and /i = 2A. Moreover, the function 
£(m + c) simplifies when evaluated at the extremal point: 



£(m + A)| 



A 

M=2A — ^ 



(3.66) 



Therefore, the expression (13.521) for the vacuum superpotential becomes: 



eff 



E 92i 

E 



#2i 

2i 



2i 



2i 



n+l i 



i v 2' 
i=1 fc=Lf J+i 



£(m + A) 



i-2fc 



n+l 



a 2, + Et E 



i=l fc=L4j-t-l 



A- 



(3.67) 



Note that by contrast to the previous example, the effective superpotential now has the 
form of a finite series. 

As we discuss in the next section, we should expect to recover this result by factor- 
izing the Seiberg-Witten curve for J\f = 2 Yang-Mills with Nf massive hypermultiplets. 
In section l3~2. ll section we will solve the factorization problem for general Nf and verify 
the equivalence of the resulting vacuum superpotential for the case Nf = N c . 

Other examples can be treated similarly by solving the extremization condition 



(13.561) to find the extremal size of the cut in the spectral curve, and substituting the 
result into (I3.52I) . These equations are exact, in that they receive no further quantum 
corrections, but in general they can only be solved as a series expansion. 



3.2 Seiberg-Witten curves and supersymmetric vacua 

In previous sections we studied the vacua of the J\f = 1 gauge theory directly. These 
results descend from the structure of the underlying Af = 2 theory one obtains by setting 
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Wtree — 0, and we turn our attention now to the J\f = 2 U(N) gauge theories with Nf 
fundamental hypermultiplets. 

As is well-known, the vacuum structure of J\f = 2 gauge theories are described by 
a fibration of a Riemann surface (the Seiberg-Witten curve) over the moduli space. At 
points in the moduli space where the curve degenerates, physical degrees of freedom 
(monopoles, dyons or W-bosons) become massless. 

For example, the Seiberg-Witten curve of Af = 2 U(N) or SU(N) pure gauge theory 
is the genus N — 1 hyperelliptic curve 

y 2 = P N (xf - 4A 2Ar (3.68) 

where Pn{x) = x N + J2iLi SiX N ~ k , with s\ = for the SU(N) curve, and A is the 
dynamically generated scale of the gauge theory. 

Written in J\f = 1 language, the effective superpotential for the J\f = 2 theory 
in the neighborhood of a point where / monopoles simultaneously become massless is 
HSW94bllSW94a1 

i 

W(M m , M m , u p , A) = J2 M m M m a Dtm (u p , A) (3.69) 

m=l 

where M m are the monopole hypermultiplets, a_o,m are the periods of the Seiberg-Witten 
curve that determine the monopole masses, and u p are the gauge-invariant curve moduli 

u = -Tr$ p (3.70) 

P 
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that parameterise the vacua of the M = 2 theory. After breaking to J\f = 1 by the 
addition of a tree-level superpotential, the Intriligator-Leigh-Seiberg linearity principle 
|ILS94| implies that the exact superpotential becomes ICIV01I 



The equation of motion for the monopole fields imposes that an,m = 0. This is true 
iff the corresponding B-cycle of the Seiberg-Witten curve degenerates, therefore the 
vacua of the gauge theory are associated to a "factorization locus" in the moduli space 
of the Seiberg-Witten curve, where I cycles of the Seiberg-Witten curve simultaneously 
pinch off to zero volume. The equation of motion for the u p then implies that there is 
a nonzero monopole condensate in the confining J\f = 1 vacua, i.e. confinement of the 
J\f = 1 theory is associated to monopole condensation. 

The maximally-confining vacua correspond to the point in the J\f = 2 moduli space 
where all iV — 1 monopoles become massless, and the Seiberg-Witten curve degenerates 
completely to genus 0. 

After evaluating (13.711) at the factorization locus, the exact effective superpotential 
then becomes 



Thus, evaluation of the effective superpotential is equivalent to solving the factorization 
of the spectral curve. Once we know the moduli (u p ) at the factorization locus we can 
immediately read off the effective superpotential corresponding to any given Wf ree using 



W(M m ,M m ,u p , A,g p ) = ^ M m M m a Dim (u p , A) + ^# p Mp ( 3 - 71 ) 



m=l 




(3.72) 



dT72h . 



The factorized Seiberg-Witten curve can be written as 



y 2 = G 2 (x)F 2 ( N -i){x) 



(3.73) 
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where the I double roots of the factorization correspond to the collapsed cycles. Since 
these collapsed cycles correspond to monopole fields that are frozen to a particular 
vacuum expectation value, they are no longer dynamical and the double roots can be 
dropped from the factorized curve, giving a "reduced curve" that describes the remain- 
ing low energy J\f = 1 dynamics. This curve is to be identified with the J\f = 1 curve 
y 2 = W'(x) 2 — fn-i(x) studied in section |3~TTT1 

For J\f = 2 U(N) or SU(N) pure gauge theory, the factorization of the curve is 
achieved as follows IDS95bB : 

y 2 = P N (x) 2 -AA 2N 

= 4A 2N (T N (x) 2 - 1) (3.74) 

where T N (x) are the Chebyshev polynomials of the first kind, defined by 

T N (x = cos(0)) = cos(iVT) 

l-l 

N y4 (~l) r 

~2 ^ N -r 

r=0 

which gives the expansion of cos(N9) in terms of cos(0). In other words, by tuning 
the parameters s k of the curve (equivalently, the gauge-invariant moduli u k = ^Tr$ fe , 
which are related to the s k via ku k + ks k + Yn=i ^ u i s k-i = 0), we can obtain Pn(x) = 
2A N T N (^), therefore 

P N (x) 2 -A 2N = A N (cos 2 (N6) - 1) = A N (sm 2 (N6)) 




N-2r 



(3.75) 
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where Un{x) are the Chebyshev polynomials of the second kind, given by 



LfJ 



Un(x) = £(-l) r 



n — r 



(2x 



,n—2r 



(3.77) 



r=0 



From (13.751) one can read off the values of the Sk in this vacuum. To convert to Uk 
we use the product form 



T N (x) = 2"- 1 l[(x - cos( l ' fc 2iv lj7r ) = ]J(x - **) (3.78) 

fc=l k=l 



with 



A' 



1=1 



Expanding the power sum for SU(N) gives 



Uk= \ 



k odd 



A k k even 



k/2 



and therefore we have the effective superpotential 



(3.79) 



(3.80) 



W 



A: 



(3.81) 



The result for Z7(iV) may be obtained from (13.811 ) by shifting x — > x + Ui/N = x — <p 
in (13.791) to account for the non-zero trace of $, where the equality follows since we are 
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in a maximally-confining U (N) vacuum, for which classically ($) = diag(0, (ft, . . . , (p). 
Explicitly, for the maximally-confining U(N) vacua, 

(;)(:h- 

If we wish, we can rewrite this expression in terms of the gaugino bilinear S, by per- 
forming a Legendre transformation with respect to the corresponding source log(A 2JV ) 
(i.e."integrating in S") IFer()3al : 



A 2N 

W(4>,g p ,S,A 2 ) = 22g p u p (<t>,A 2 = y) + Slog' 



yN 

p>i y 



p \ ( 2q 



LP/2J 

p >i P q =o \ 2q I \ q 



A 2N 

+S\og(—) (3.83) 

yN 

Of course, this form of W does not contain any additional information, but it is useful 
for comparison with the other techniques we discuss. For example, when we study the 
relationship between effective superpotentials and integrable systems in chapter 0] we 
will recover this expression from the Lax matrix of the affine Toda system. 

3.2.1 Factorization of the Seiberg-Witten curve for Nf > 

The Seiberg-Witten curve for Af = 2 gauge theory with < Nf < 2N C fundamental 
hypermultiplets is ||DKP97 1 

y 2 = P Nc (x) 2 - 4A 2N <~ N f l[(x + Tru) (3.84) 

i=i 
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where m ; are the bare hypermultiplet masses. When Nf > N c there is an ambiguity in 
the curve, and a polynomial of order Nf — N c in x (multiplied by appropriate powers of 
A to have well-defined scaling dimension n) may be added to Pn c (x) without changing 
the J\f = 2 prepotential. For comparison to the results of section U. 1.61 we will mainly 
be interested in the case Nf = N c for which the ambiguity in Pn c ( x ) appears at constant 
order and is proportional to A^. 

The curve (13.841) can be scaled to recover the Nf = curve (13.741) by taking the limit 

A^O, m^oo, A 2N <~ N f Y[ mi = A 2N c (3.85) 

with A finite. Note that the latter identification is the scale-matching relation of the 
theories above and below the mass scale of the fundamentals. 

We now show how the factorization using Chebyshev polynomials can be general- 
ized to the hypermultiplet curve (13.841) (this problem has been studied indirectly using 
matrix models in [DJ03b|). Define the functions 

PnM = J>cos((^c-*")0) 

i=0 
i=0 

Then 

Pk-Qk = ^^ + 2^^(cos(^)cos(^)+sin(^)sin(^)) 

i i^j 

= J2^ + 2 J2 - 3)0) = R Nf (0) (3.87) 

i i^j 
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Therefore the equation 

Pn c ~ R-N f = Q% c (3.88) 

gives the desired factorization of the Seiberg-Witten curve by setting cos(#) = ^ for 
U(N), or cos(#) = for SU(N), where the shift is needed to cancel the x N ~ l term 
in Pjv(x). The parameters I/; are related to the fundamental masses although the 
relations are polynomial in general. 

This expression simplifies dramatically when Nj = N c , mi = to, and we find 



N 



N 

i 



P N ~ l cos(i9) 




(3.89) 



with P = A/A, where A is the scale of the theory with flavors, and A 2 = mA is the 
parameter defined above that corresponds to the dynamical scale of the theory in the 
limit where the fundamentals have been scaled out completely. If we choose a limit 
where the fundamental masses become very large compared to the scale A, i.e. such that 
P becomes a small parameter, then the curve can be treated as a small deformation of 
the Nf = curve. 

After some algebra, we obtain the following expression for Pn(x): 



W-UX + fi { N ) ('-") ^ (3,0, 

i=l \ % J r=0 \ T J 



78 



where A = for U(N) and A = A for SU(N) to cancel the first subleading power of 
x. We can resum this expression to extract the s fc [Mor56|. For U(N) we find 



sn 



-i = A""' £ (j + 2r) 



AT 



r=0 



2r / -? + r 



A 



2r 



+ 2A%- (3.91) 



and for SU(N) we find 



Sat 



-i = a"-'E 



i=i 



+2A"S ji 



E 

r=0 



i — r 



*-2r \ ( A 



2r 



(3.92) 



We now compare to the results obtained in section l3~ 1 .4l based on period integrals of 
the J\f = 1 curve. Recall that for Nf = N c we obtained the expression (13.671) 



W eff =N c 



i=l 



Ef I" |A 2, +Ef E 



2i 



n+l 



i=l fc=LfJ+l 



(3.93) 



From this expression can be read off the values of the gauge-invariant moduli (u k ) = 
Note that our result has the form of a finite series expansion in (3, and in the limit 



(3 = we recover the superpotential (13.461) of the Nf = theory. The u k are related to 
the curve parameters s k via the Newton formula 3 



fc-i 



ku k + ks k + 2j ifHSk-i = 



(3.94) 



i=i 



3 This footnote inserted to see if anyone notices it. 
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As in section l3~2l the SU(N) moduli Uk may be obtained from the U(N) by shifting 
away the trace: 

N 

fi * = E& - (3 - 95) 

»=i 

Expanding the powers in ( 13.951 ) one finds 

** = i I E(^)*-^ ( * j ^ + (3-96) 

We have verified in a number of cases that the Uk associated to the Sk (13.911) obtained 
from the factorized Seiberg-Witten curve agree with the values calculated from the 
superpotential ( 13.931) . up to a physically irrelevant sign A — > —A (which can be absorbed 
into the conventions used to define the Seiberg-Witten curve (13.841) ) and the ambiguity 
in the top modulus at order A^. 

For example, the factorization for the first few U(N) curves is achieved by: 

U(2): 

p 2 {x) = x 2 + 2xA - 2 A 2 + 2 A 2 

ui = -2A, 

u 2 = 2A 2 (3.97) 
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t/(3): 



P 3 (x) = x 3 + 3x 2 A + x(-3A 2 + 3A 2 ) - 6A 2 A + 2A 3 

ui = — 3A, 

u 2 = 3(A 2 + iA 2 ), 

^ = 3(-A 2 A-w\ 3 ) (3.98) 



1/(4): 



p 4 ( x ) = x 4 + 4x 3 A + x 2 (-4A 2 + 6A 2 )+x(-12A 2 A + 4A 3 ) 
+2A 4 - 12A 2 A 2 + 2A 4 
mi = -4A, 
u 2 = 4(A 2 + iA 2 ), 
m 3 = 4(-A 2 A 2 -iA 3 ), 

w 4 = 4(^A 4 + A 2 A 2 ) (3.99) 

which can be compared to the Uk read off from (13.931 ): 

mi = iVA, u 2 = N(A 2 + ±A 2 ), us = N(A 2 A + ±A 3 ), 
M 4 = A r (|A 4 + A 2 A 2 + |A 4 ) 

3.3 SOISp gauge groups 

In chapter |5] we will describe our work llACH + 03l on using matrix models to compute 
effective superpotentials for SO and Sp gauge theories. Here we give a short discus- 
sion of field-theoretical aspects involving the Seiberg-Witten curve. The Seiberg-Witten 
curves for J\f = 2 pure Yang-Mills theory with SOISp gauge groups were found by 
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IIDS95a[lBL95[lAS96ll . For a rank-r gauge theory, the spectral curve is a genus r hyper- 
elliptic curve, of the form 

y 2 = P2r + 2{x i {<P l }) 1 (3.100) 

where P 2 r+2 is a polynomial of degree 2r + 2 in the x that also depends on the moduli 

The SO and Sp spectral curves can also be written as a genus 2r — 1 curve, 

y 2 = P 2r {x\ {&}), (3.101) 

which is therefore symmetric under the Z 2 action x i— > —x and is a double cover of the 
genus N curve (13.1001) via this map. This is the form of the curve we will work with. 
Because of the Z 2 symmetry, each point is paired with its image; this will be important 
when we discuss the matrix models for SO and Sp gauge theories in section 1531 since 
the matrix model eigenvalues live on the factorization of this curve, and therefore also 
come in pairs. 

J\f = 2 supersymmetry may again be broken to J\f = 1 by an appropriate 
gauge-invariant superpotential term for $. Because the trace of odd powers of matri- 
ces in the Lie algebra of SO(N)/Sp(N) vanishes, the superpotential deformation for 
SO(N)/Sp(N) only includes polynomial terms of even degree: 

n+l 

w tree m = Y,^n^ k )- (3.io2) 

k=l 

A superpotential Wf ree of order 2n + 2 breaks the gauge symmetry down to a direct 
product of n + 1 subgroups, e.g.: 

SO(N) -> SO(N ) x U(Nx) x ... x U (N n ) , (3.103) 
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where N = N + 2N X + ■ ■ • + 2N n . 

We saw in section l3~2"l that the supersymmetric vacua of the A/" = 1 theory require r — 
n mutually local monopoles to simultaneously become massless and condense, leading 
to confinement of the gauge theory. Imposing this condition is therefore equivalent to 
the factorization HCV02I 



where pi ^ pj, q { ^ qj for i ^ j. On this locus we then obtain (after discarding the 
terms corresponding to the non-dynamical condensed monopoles) the reduced spectral 
curve 

2n 

y i = Y[(x 2 -q*), (3.105) 

i=i 

which has genus 2n— 1. This curve parameterizes the J\f = 2 vacua that are not lifted by 
the deformation to J\f = 1 (I3.102I) . Notice that the curve is still invariant under x \—> —x: 
this implies that the branch points come in pairs: (— (&). This reduced spectral curve 
is identified with the curve (13.101 ) arising from string theory discussed in section |3~T1 

The factorization problem was solved by JTO03| along the lines of the SU(N) dis- 
cussion in section l3~2l however as discussed in section l3~TTT1 the effective superpotential 
of the J\f = 1 gauge theory can also be obtained from the periods of (13.1051) . It will take 
the form 



r—n In 



y 2 = W^-pl) 2 \{^-il 



(3.104) 



t=l j=l 




(3.106) 
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where AmSi are the periods of the meromorphic 1-form y dx around the A-cycles of the 
spectral curve, flj the corresponding periods around the B-cycles, and Ni is 



N 



N SU{Ni), 
1 SO(Ni), 



(3.107) 



k N, + 1 Sp(N t ] 



In chapter|5]we will see how the shift Ni i— > N emerges from a subleading correction 
to the effective superpotential, obtained using matrix integral techniques. 
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Chapter 4 

Integrable Systems and Af = 1 Vacua 

In this chapter we investigate the relationship between J\f = 1 superpotentials and inte- 
grable systems. Integrable systems are known to underly the low energy structure of 
gauge theories with J\f = 2 supersymmetry (|MW96|, see |DP99| for a review), and 
this underlying structure again survives soft supersymmetry breaking to TV = 1 to gov- 
ern the effective superpotential. The origin of these integrable structures in gauge theory 
is still incompletely understood. In this chapter we obtain some new details about the 
correspondence, and it is an open problem to understand the results in a more general 
context. 

In section 13.2.11 we obtained simple combinatorial formulae for the moduli Uk of 
the M = 2 Seiberg-Witten curve at the maximal factorization locus. For Nf = 0, 
these combinatorial formulae are encoded in the traces of powers of a particular matrix, 
namely the scalar component of the adjoint field $, evaluated in the vacuum of interest: 

(u k ) = ^Tr$ fe (4.1) 

The connection to integrable systems is via this matrix $, which is identified with 
a Lax matrix for the associated integrable system. The Lax matrix completely defines 
the dynamics of the integrable system. The integrable system associated to pure J\f = 2 
Yang-Mills is the periodic Toda chain [ MW96 1 , and the known Lax matrix of this system 
can indeed be identified with ($). An algorithm was presented by [BdBDW04| for 
computing the effective superpotential of Af = 1 gauge theory with an adjoint chiral 
superfield, using this Lax matrix. 
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The integrable system associated to J\f = 2 SQCD (i.e. Nf ^ 0) was uncovered in 
[GMMM96a 1 GMMM96b|, and is a particular spin chain system. However the known 
Lax pair of this system is written in transfer matrix form as a chain of 2 x 2 matrices, 
for which the connection to J\f = 1 superpotentials is less direct since this form of the 
Lax matrix does not have an obvious physical meaning 1 . Therefore, it would be useful 
to find another Lax pair for this system that takes the form of a single matrix, similar to 
the Nf = case 2 . We studied this problem in [KW03|, and in section l4~2l we present 
the N c x N c matrix (<&) that encodes the (uk) in the maximally-confining vacua, which 
is identified with a particular equilibrium value of the Lax matrix for the associated spin 
chain. 

We begin this section by proving that the superpotential calculation of [BdBDW04| 
using the integrable structure of the J\f = 2 gauge theory, yields the same result in the 
maximally-confining phase as (I3.46I) . (13.831) obtained from the Nf = period integral 
and factorization calculations [K W03I . 

'See however the recent work [HO04|. 

2 Note that a given integrable system may have more than one Lax pair, and the matrices may even be 
of different rank, so we should not be discouraged from looking for a new Lax formulation. 
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4.1 The periodic Toda chain and N = 1, Nf = vacua 

The integrable system associated to pure J\f = 2 Yang-Mills theory is the periodic Toda 
chain, which has Lax matrix: 



/ 



1 






yi o 

02 2/2 







z 




Vn-i 

4>N 



\ 



(4.2) 



y y N / z ... 1 07V J 
where i5 ?/j are the dynamical position and momentum variables of the integrable sys- 
tem, whose precise definition will not be important for us (see [DP99| for a review), and 
z is a "spectral parameter", an auxilliary variable not associated to the physical system. 
The conserved quantities (Hamiltonians) of the Toda system XJ k = ^TrL fc are associated 
to the gauge-invariant polynomials w fc = ^Tr$ fe that parametrize the moduli space of 
the M = 2 gauge theory. The spectral curve of the Lax system is defined by 



, 2N -1^ 



det(x.J - L) = Pn(x) + (—l) (z + A z 







(4.3) 



where Pn are the polynomials defined in section l3~2l Under the change of coordinates 



y = 2z + (-l) N P N (x 



(4.4) 



the spectral curve becomes 



Pn(x) 



4A 



2JV 



(4.5) 
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which is the standard form of the Seiberg-Witten curve of J\f = 2 U(N) Yang-Mills 
theory. Therefore, when we deform the J\f = 2 theory by turning on a tree-level super- 
potential 

n+l 

" (4.6) 



n+l 

W = 9iUi 

i=l 



the analogous quantity in the Toda system is the corresponding function of the conserved 
quantities U^. The essence of the proposal of [BdBDW04] is that evaluating W(L) gives 
the exact effective superpotential of the theory 3 . The factorization of the spectral curve 
at the points corresponding to J\f = 1 supersymmetric vacua translates in the integrable 
system to equilibrium configurations that are stationary under the Hamiltonian flows 
generated by the U k [Hol03|. 

We will now obtain the explicit form of Wi ax for a given Wf ree and recover the result 
in section I3T21 For this purpose the form of the Lax matrix (14.21) is slightly awkward to 
work with, because the z entries are not on the same footing as the other variables. To 
rectify this, conjugate L by diag(l, z 1 ' 1 *, z 2 l N , . . . , z N ~ l l N ) to bring it into the form: 



/ 



A/N 






7TJN 







y 

7TW 















1/2V \ 








y 

.l/N 



(4.7) 



3 When the superpotential Wt ree contains terms of degree N or higher, the spectral parameter z that 
appears in the Lax matrix J4.2i will appear in the Uk- However, in the quantum Af = 1 gauge theory 
these moduli are ambiguous because the operators Tr$ fc , k > N receive quantum corrections, and the 
resolution proposed in |BdBDW04| was that all occurrences of z in the Lax superpotential W(L) should 
be discarded at the end of the computation (alternatively they can be suppressed to arbitrarily high orders 
by embedding U(N) C U(tN)). 
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where S is the N x N shift matrix, satisfying S = I. 
Therefore, 




(4.8) 



where in the second line we have used the fact that the terms can only appear on the 
diagonal if 2m — I = Na, a G Z. Suppressing powers of z whenever they appear, we 
obtain 

W lax = N £ f £ ( " ) f 2 " ) r- 2 V + S log(^) (4.9) 
p>i F q =o \2q J \ q J » 

which recovers the expressions (13.46b . (13.831) obtained using exact field theory tech- 
niques, and by evaluating period integrals. 

4.2 Results on a new Lax matrix for Nj = Af c 

The connection between J\f = 2 gauge theories and integrable systems can be summa- 
rized by identifying the matrix-valued field $ of the quantum gauge theory with a Lax 
matrix for the integrable system. Therefore, if we can evaluate ($) in a given vacuum, 
we know the value of the Lax matrix in an equilibrium configuration of the integrable 
system. Knowing the values of the moduli in the particular J\f = 2 vacuum gives 
iV c equations for the matrix ($), which is enough in principle to determine (<&) up to 
gauge transformations. 
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In the previous section, we showed how evaluating the Toda Lax matrix in a particu- 
lar equilibrium configuration (all position and momentum variables equal, i.e. (pi = (ft, 
y,i = y = A 2 ) allows us to recover the (u k ) of the factorized Seiberg-Witten curve. 
Conversely, given the we can reconstruct the Lax matrix of the periodic Toda 
chain: the (uk) may be obtained from the single matrix 4 



/ 



A 2 

A 2 

. . . 1 

... 






2A 2 



(4.10) 
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One can explicitly see from this expression how the classical value of $ = 
diag(0 1; . . . , <p N ) is deformed by quantum effects, specifically the interaction with 
the background magnetic field of the condensed monopoles, which generates the off- 
diagonal terms (this can most easily be derived via compactification to 3 dimensions, 
where the four-dimensional monopoles reduce to 3-dimensional instantons [dBH097 1). 

We follow the same philosophy for the Nf = N c vacua studied in section 13.2.11 
and identify the matrix <£> from which the expectation values of the moduli (uk) in the 
maximally-confining vacua may again be obtained by taking the trace of powers (recall 
that in this case the moduli took the form of a finite series). We therefore have a candi- 
date for a Lax matrix of the associated integrable system, which in these examples are 
spin chains ||GMMM9"6al IGMMM96b1 . 

4 The entry with coefficient 2 exists because (14. lOi does not contain the spectral parameter z (which 
does not have a physical meaning in the gauge theory), so we can absorb the entry A 2 /z of ( 14. 2\ into this 
entry. 
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We find for SU(N C ), Nf = N c and all quark masses equal, that the moduli (uk) of 
the maximally confining vacuum may be obtained from the matrix 



<<&) 



V 






A 2 


AA 2 


A 2 A 2 




A^- 2 A 2 


NA^K 2 


1 





A 2 


AA 2 




A^- 3 A 2 


(N - 1)A N ~ 2 A 2 





1 





A 2 




A^ 4 A 2 


(N - 2)A N - 3 A 2 











1 





A 2 


3AA 2 














1 





2A 2 

















1 






(4.11) 



Note that this reduces to the Toda Lax matrix (14.101) in the appropriate scaling limit A — > 
(here = for the SU(N) vacua to ensure tracelessness). It remains an open problem 
to generalize this matrix to a general vacuum and to better understand the relationship 
with the degrees of freedom of the spin chain system. 
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Chapter 5 

The Combinatorial Structure of 
Supersymmetric Vacua 

In this chapter, we first introduce matrix integrals and review how they may be solved 
in a genus expansion. The solution produces a "spectral curve", which is isomorphic 
to the Af = 1 curve (13.101) when the potential of the matrix model is identified with 
the tree-level superpotential for the adjoint chiral superfield W(&). The matrix integral 
is a generating function for Feynman diagrams that are in one-to-one correspondence 
with planar $ diagrams for the SU(N) gauge theory, although the diagrams of the four- 
dimensional gauge theory carry four-dimensional momenta in the loops whereas the 
matrix diagrams are zero-dimensional and carry no internal momenta. Nonetheless, as 
we discussed in chapter |3] this spectral curve produces the exact effective superpoten- 
tial of the four-dimensional gauge theory in terms of integrals of the resolvent along 
contours of the curve. The insight afforded by the matrix model is that it provides a 
remarkably simple perturbative expansion of this effective superpotential. 

The relationship between matrix integrals and gauge theory superpotentials was first 
observed for Af = 1 SU(N) gauge theories with adjoint matter [DV02a| and for the 
Af = 1* deformation of Af = 4 SU(N) SYM lDV()2cl IDHPKS02L IDHKS021 The 
conjecture was subsequently extended to a number of other cases including HDV02bl 
ICM03I IDHK02I IFer()3bl lBer03l IFO021 |Gop03| IDJ()3al IACFH()3bl IMcG()31 ISuz()31 
IBR031 lGor03l INSW()3al ITac03l IDNV021 IKMT03I IDST021 IFen()21 IFH031 IACFH()3al 
NSW03b|. 
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The underlying quantum field theoretical reason for the correspondence is that after 
summing the four-dimensional Feynman diagram contributions of the field $ to the 
effective superpotential, all dependence on internal loop momenta cancels and one is 
left with only the zero-momentum planar $ diagrams of the theory llDGL + 03 1 IAI VW031 
(and for SOISp gauge groups, and theories with fundamental matter, the leading non- 
planar), which are generated by the associated matrix integral. In section ISTTI we review 
the technology of zero-dimensional matrix models, the computation of the matrix model 
free energy and the gauge theory effective superpotential. 

In order to understand the combinatorial origin of the divergent period integral 
appearing in the matrix model calculation, we extend the matrix integral to include 
M x 1 and lxM vectors (corresponding to gauge theory with matter in the fundamen- 
tal and antifundamental representations). The generating function of planar diagrams 
with 1 boundary recovers the previous expression (13.481) for the quark contributions to 
the superpotential in terms of period integrals of the spectral curve. As in the geomet- 
rical analysis of chapter |3l adding Nf = 2N C vectors to the matrix potential causes the 
divergences of the period integral to cancel, and we see explicitly the role of the "quark" 
vectors in regularizing the matrix integral computation of the effective superpotential. 
Thus, the cancellation of this divergence is understood perturbatively in the gauge the- 
ory as coming from the contribution to the effective superpotential of planar Feynman 
diagrams with disk topology, in the limit where the quarks that propagate around the 
boundary of the diagram are very heavy. This result was published in IK W03I . 

Finally, we extend the analysis to gauge theories with SO and Sp gauge groups, 
which was first published in llACH + 03l . This amounts to including in the counting the 
non-orientable Ribbon diagrams with MP 2 topology. This requires an adaptation of the 
technique of higher-genus loop equations [ACKM93, Ake96]. We show that the MP 2 
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contribution to the resolvent - and hence to the matrix model free energy - has a simple 
form and is related to the genus-0 result. 

5.1 Matrix integrals and zero-dimensional matrix mod- 
els 

Consider the matrix integral 

Z = Z J dMe-^ wm (5.1) 

where M G Q is an M x M matrix, VT(M) is a polynomial, and Zq is a normalization 
factor. The integral can be rewritten in terms of the eigenvalues of M as 

. M 

Z = I fldWaWe-^™ (5.2) 

where J is a suitable Jacobian for the change of variables, and Z is fixed by the normal- 
ization of (15.2b . For Hermitian matrices, the change to the (diagonal) eigenvalue basis 
involves conjugation by unitary matrices 

M -> U f DU (5.3) 

and the change of basis produces an integral over the Haar measure on U(M), which 
gives the volume of U(M) / U(1) M and fixes Z . The Jacobian is given by J = 
A 2 ({Aj}), where A is the Van der Monde determinant 

A({A,}) = - Aj) (5.4) 

i<j 
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The action for the matrix eigenvalues is then 

1 M M 

s{\) = — w ^ + E E lo ^ - A *) (5 - 5) 

^ s i=l i=l «<j 

We will mostly be interested in the 't Hooft limit 

M -> oo, 5 = = con*?. (5.6) 

Recall from section I3TT1 that this is also the limit in which gravity decouples from the 
string theory. We wish to evaluate the matrix integral in this limit as a perturbative 
expansion around a saddle point (classical vacuum). Such a classical vacuum is given 
by a distribution of the eigenvalues of M among the critical points {xi] of the function 
W(x). We denote the number of eigenvalues at the critical point Xi by Mj and define 
the corresponding 't Hooft couplings Si = g s Mi. 
The free energy T of the matrix model is given by 

T = log Z = TP ert + T non -P ert (5.7) 

where pP ert - comes from evaluating the integral perturbatively (as usual, it is given by 
the sum of connected Feynman diagrams), and p non ~P ert - i s a non-perturbative contri- 
bution that will be determined later. 

Consider the propagator for the Hermitian matrix model: it has a group theoretical 
factor 

(MyM«) ex 5 u 5 jk . (5.8) 

The propagator and interaction vertices may be represented in double-line notation, see 
Figure l5~T1 We can now expand the free energy perturbatively around a given classical 



95 



a) 



b) 



Figure 5.1: Feynman rules for the Hermitian matrix model: a) propagator, and 
b) sample quartic vertex, giving the perturbative expansion in terms of "ribbon 
graphs". 

vacuum in terms of ribbon graphs, where the edges of the ribbon correspond to eigen- 
values Aj of the matrix M. Thus, each closed loop of a ribbon graph edge contributes 
a factor of M, = Si/g s , the number of eigenvalues on the 2 th critical point. From the 
overall normalization of the action (15.51) . it is clear that each vertex of the diagram con- 
tributes a factor of l/g s and each propagator contributes g s . Thus the overall power of 
9s is 

9 P s~ v ' 1 = 9S X = 9 2 r 2 (5-9) 

where p is the number of propagators, v the number of vertices, / the number of index 
loops and \ — 2 — 2g is the Euler characteristic of the Riemann surface with minimal 
genus g on which the diagram may be drawn. Therefore the perturbative free energy has 
a topological expansion 

F(gs,Si) =F n P(S l )+Y,9s 2g - 2 J r g (S l ) (5.10) 

9 

In the 't Hooft limit (15.61) the planar (sphere-topology) diagrams dominate 1 . 
'Higher-genus contributions correspond to gravitational corrections from string theory. 
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Note that we are distinguishing the rank M of the matrix from the rank N c of the 
four dimensional gauge theory we will soon make contact with. As we have mentioned, 
it is only the planar diagrams of the four-dimensional gauge theory that contribute to the 
effective superpotential even for finite N c , so the large M limit taken in the matrix model 
is an auxilliary step designed to isolate the planar diagram contributions to the value of 
the matrix integral. The explicit dependence on M is hidden by rewriting M = S/g s , 
and S is identified with the gaugino bilinear superfield of the gauge theory. 

A non-perturbative contribution to the free energy comes from the residual gauge 
invariance that exists when two or more eigenvalues coincide. When Mj eigenvalues are 
distributed in the 2 th critical point of the potential, the matrix integral is invariant under an 
additional niLi U(Mi) gauge symmetry. Thus, the path integral includes the orbit of the 
solution under this group, so the free energy of the matrix integral receives an additional 
contribution from the logarithm of the volume of these gauge factors [Mor95 , OV02]: 



This point was unclear in much of the literature, and the volume contribution was often 
confused with the normalization Z of the matrix integral (15 .lb . However, this would not 
give the correct contribution in vacua with broken gauge symmetry, and the logarithm 
contributes with the opposite sign relative to the perturbative terms, which can be ruled 
out by an explicit evaluation of the free energy using the techniques we discuss below. 

The asymptotic expansion of the volumes are worked out in Appendix lAl For U(M) 
we obtain 




(5.11) 
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Changing variables using M = S/g s and extracting the leading term in g s , we find that 
the leading-order non-perturbative contribution to the free energy is 



o 



np 



S -\ og {S/g s )+ l -S 2 \og2K + - A S 2 



(5.13) 



These non-perturbative matrix model contributions coincide with terms that describe 
non-perturbative physics in the gauge theory, namely the Veneziano-Yankielowicz 
superpotential, which is associated to the strongly-coupled gauge dynamics. We will 
discuss this more later. In principle there may be other non-perturbative contributions to 
the free energy. In the limit of large M we can evaluate the matrix integral by the method 
of steepest descent, which will allow us to compute the free energy directly. It can then 
be verified in examples that the leading contribution to the free energy reproduces the 
perturbative and non-perturbative terms discussed above. 

We will now discuss the solution of the matrix model in the eigenvalue basis. From 
the eigenvalue action (I5.5I) . the equation of motion for a single eigenvalue Aj is 



After multiplying by 1/ (Aj — z) and summing over i, equation (15.141 ) becomes 




(5.14) 



It can be solved by introducing the resolvent 




(5.15) 



u 2 (z) - g s u'{z) - W'{z)u{z) - -J{z) = 



(5.16) 



where 




(5.17) 



98 



Equation (15.161) is the classical loop equation. In the 't Hooft limit (I5.6I) . the second 
term in (15.161 ) can be neglected, and performing the change of variables 



y(z) = 2u(z) - W\z) 



(5.18) 



it reduces to 



y\z) = w\zf + f(z) = o 



(5.19) 



The coefficients in f(z) are as yet undetermined. Note that this curve has the same 
form as the curves discussed in the previous chapters, if the matrix potential iy(M) 
is identified with the tree-level superpotential for the adjoint chiral superfield The 
recovery of this curve from the matrix model is a signal that the matrix model is related 
to the four-dimensional gauge theory, since the same curve also emerged from string 
theory when we studied geometric engineering of the gauge theory in section I3TTI 



Using (15.181) and (15.191) the equation for the resolvent yields a formal solu- 
tion |DFGZJ95| 



where the branch of the square root is fixed by the requirement that the resolvent have 
asymptotic falloff uj ~ S/x, which vanishes in the classical limit S — > 0. The resolvent 
is thus expressed in terms of the n unknown coefficients that appear in the polynomial 
f(x) defined in ( I5.17I ). From the form of the solution, it is clear that the resolvent 
has square root branch cuts around the critical points of the matrix potential iy(M). 
Physically, the eigenvalues sitting at the critical points feel a Coulomb repulsion from 
the logarithmic term in the eigenvalue action (I5.5I) . and spread out from their classical 
values to form the cuts. 




(5.20) 
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In the large M limit the distribution of matrix eigenvalues 



p(A) = ^5(A-A l ) 



(5.21) 



becomes continuous. In terms of p(\) the resolvent can be rewritten as 




(5.22) 



which implies that 



P(A) 



1 



(u(\ + iO) -w(A-iO)) 



1 



(y(\ + iO)-y(\-iO)). (5.23) 



2nig s 



s 



i.e. the eigenvalue density is given by the discontinuity of the resolvent across its branch 
cuts. The 't Hooft parameters associated to the number of eigenvalues in the i th cut are 
then given by 



The function y(x) contains the singular part of the resolvent. It can also be written 

as 



where S is the action, the derivative of which gives the force acting on an eigenvalue. 
Now, if the number of eigenvalues on the 2 th cut is varied by taking an eigenvalue to 
infinity along the non compact Bi contour of the Riemann surface (15.481) . the change in 
the free energy JF of the matrix model is given by the line integral of the force along 
this contour: 




(5.24) 



(5.25) 




(5.26) 
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This is a differential equation that determines JF , the leading (genus 0) contribution to 
the free energy of the matrix model. By equation (15.251 ) it is apparent that (15.261 ) gives 
the action for an eigenvalue to tunnel from the cut to infinity. 

Note that the contour integral (15.261) is again logarithmically divergent, as we saw 
when we encountered the same integral in the context of factorized Seiberg-Witten 
curves. As we will see in the following section, it can again be understood in the matrix 
model by introducing 2N C vectors into the matrix potential, whose planar combinatorics 
cut off the integration contour and render the integral finite. Thus, the matrix model 
provides a simple intuitive interpretation of the spectral curve and the related period 
integrals, in terms of the dynamics of matrix eigenvalues. 

The fact that we have recovered the same curve from the matrix model suggests 
that the matrix model is related to the string theory (recall that the curve was obtained 
from the Calabi-Yau compactification manifold) and to the four-dimensional gauge the- 
ory that it engineers. Indeed, this is the case, as first discussed in the seminal work 
IDV0 2al: the action of B-type topological strings on the Calabi-Yau spaces of section 
13-ll reduces to the matrix models, and at the same time computes the gauge theory effec- 
tive superpotentials. 

Since the spectral curve (15.191) and meromorphic 1-form y dx are the same as those 
obtained from the Calabi-Yau geometry discussed in section l3Tl the genus free energy 
Tq of the matrix model is identified with the prepotential of the Calabi-Yau. In other 
words, the large M solution of the Hermitian 1 -matrix model (15.11) computes the pre- 
potential of Type IIB string theory on the associated Calabi-Yau manifold 2 . Therefore 
as discussed in section I3TTI the superpotential of the gauge theory is given in the matrix 

2 This is also true for other types of matrix integral, although there are relatively few matrix integrals 
that can be solved exactly. 
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model by (E3 



(5.27) 



Using the result (15.131) for the non-perturbative contribution to the matrix model free 
energy we find 



where the additional logarithms in (15.131) have been absorbed into the definition of the 
cutoff 3 . Thus, the residual gauge symmetry of the matrix model vacua precisely gives 
rise to the Veneziano-Yankielowicz superpotential, which is associated to the strong- 
coupling gauge dynamics of the four-dimensional gauge theory. 

Moreover, the planar diagrams of the matrix model, which contribute to the per- 
turbative expansion of the free energy, are in 1-1 correspondence with planar Feynman 
diagrams of the gauge theory. After summing these gauge theory Feynman diagrams, all 
dependence on four-dimensional loop momenta cancels, and their contribution is effec- 
tively zero-dimensional IDG L + 03l . Thus, the matrix model can be used to compute the 
effective superpotential of the gauge theory. 

In section 13.1.21 we already evaluated the contour integral (15.261) for the simplest 
(Gaussian) matrix model (i.e. gauge theory tree level superpotential W($) = |mTr$ 2 ), 
and recovered the Veneziano-Yankielowicz superpotential together with correction 
terms that are understood as coming from the regulator superfields. Equation (15.261 ) 
can be evaluated for arbitrary 1-cut matrix models (all eigenvalues in one critical point), 
using the techniques presented in chapter |3] For multi-cut matrix models the calcula- 
tions become more difficult, but still tractable in some cases. 

3 Strictly speaking we have not yet argued for the need to introduce a cutoff into the matrix integral; 
we will justify this point later. 




(5.28) 
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As we have discussed above, the matrix model provides a simple alternative to evalu- 
ating the contour integrals: for a given vacuum distribution of eigenvalues we can simply 
enumerate planar Feynman diagrams of the matrix model to obtain the perturbative free 
energy to the desired order, and add to it the non-perturbative volume contribution from 
the residual gauge symmetries of the vacuum. Summing the perturbative expansion to 
all orders is equivalent to evaluating the integral (15.261) . 

Comparing to the four-dimensional J\f = 1 gauge theory, we see explicitly the 
origins of the Veneziano-Yankielowicz superpotential for the strongly coupled gauge 
dynamics, as well as the perturbative corrections from planar $ diagrams of the gauge 
theory. This provides an elegant physical insight into the nature of the effective super- 
potentials calculated using the techniques of chapters |3] and H which involve many of 
the same calculations, but whose physical origins are less clear. 

5.2 Matrix models for adjoint and fundamental matter 

The matrix model discussed in the previous section corresponds to J\f = 1 gauge theory 
with an adjoint chiral superfield. As in the geometrical analysis of section 13.1.11 we 
encountered a logarithmic divergence in the contribution to the effective superpotential, 
which needed to be regulated by introducing a cutoff. In the previous discussion we 
understood this cutoff as coming from 2N C fundamental chiral superfields, which sub- 
tract the divergence at infinity of the integral and replace it by a cutoff equal to the mass 
of the fields. 

The same analysis can be carried out in matrix language. In other words, consider 
the matrix model with potential 

W tree = W M {M) + J2 thQiQ* + 9iQMQ i (5-29) 

i=l 
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where M is an M x M Hermitian matrix, Qj are 1 x M vectors and Qj are transpose M x 
1 vectors. The "Yukawa couplings" can again be set to g = 1 by a rescaling of Q and 
Q. This theory has been studied for W M (M) = ±mTrM 2 in lACFH03bllBIN+03l . We 
will derive the solution to the matrix model for a general Wm using the combinatorics 
of planar diagrams, focusing on the contributions of the vectors Q to the free energy. 

As before, the matrix integral has a topological expansion, and the contributions at 
large-M now come from planar diagrams with and 1 quark boundary: 



g,h 

where g is the genus and b the number of boundaries, and we again recognise the Euler 
characteristic \ — 2 — 2g — b. Extending the result from the previous section, the 
superpotential is given by IIDV02cl lA€FH03bl 



Contributions to the first term come only from $ self-interactions, so their combinatorics 
are the same as for the theory without quarks. Diagrams with one external boundary can 
be counted by decomposing the counting problem into two parts: the combinatorics of 
the $ diagrams on the interior of the disc, and the combinatorics of the boundary of the 
disc. 

The first problem is equivalent to counting the planar n-point Green's functions 
G n (gi) of the theory without quarks (i.e. planar $ diagrams - possibly disconnected 
- with n external $ legs). This problem was solved in [BIPZ78|, as follows: 

By definition, 




(5.30) 



W(S) = N c 



(5.31) 



dS 




(5.32) 
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where the second equality follows from the change of variables from the matrix integral 
to the eigenvalue basis and a, b are the endpoints of the eigenvalue branch cut. In other 
words, the sum of the planar Greens functions at each order are given by the correspond- 
ing moment ofy(X). 

The generating function for the Greens functions is 

oo 

4>{j) = Y,3 k Gk = -w{-) (5-33) 
1 1 

k=0 J J 

where the second equality is given in terms of the resolvent 

"(A) = \{W'{\) - VW'(\y + f n -i(\)) (5.34) 

by summing the geometric series in A coming from (15.32ft . and converting the integral 
to a contour integral. We also use the previous results that the eigenvalue density p(A) is 
equal to the discontinuity in j^uj(X) across the branch cut (see (15.231) ). and has asymp- 
totic behavior uj(x) ~ S/x as x — > oo. 

To include the combinatorics of the boundary requires multiplying by ( ' fc fc , 1 ' >! = \ at 
order k in the expansion of G, to take into account the (k — 1)! distinct ways to connect 
a boundary quark with a leg of the internal Greens function 4 , and the i coming from 
the expansion of e s to order k. The factor | can be incorporated into ( 15.331) simply by 
integrating it: 

4 At first sight, it looks like an arbitrary connection of a boundary leg to an internal leg can make 
the overall graph non-planar, however we can always perform a corresponding crossing operation on the 
internal part of the diagram to undo this non-planarity. 
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n (i) = /^(j)*' 

= ™ y (W(^) - VW'ix) 2 + f n - 1 (x))dx (5.35) 

where we have changed variables x = j and used the definition of 

The factors of j count the number of external legs of the Greens function; therefore 
terms of order j k are associated to k powers of the Yukawa coupling g, and k quark 
propagators jj to connect up the k external quarks on the boundary. Therefore the one- 
boundary contribution to the matrix integral free energy is given by 



1 f°° 

Ai = -o / (W'(x) - y/W(x) 2 + fn-i(x))dx 

poo 

= — I uj(x)dx (5.36) 

and the contribution to the effective superpotential (15.311) of the planar diagrams with 
1 boundary precisely recovers the previously claimed result ( 13.481) . In the same way as 
in section U. 1.41 when there are 2M vectors the integral ( 15.361) combines with (15.261) to 
cancel the divergence at infinity, and cut off the integral at the value of the vector masses. 



5.3 Matrix Models for SO/Sp Gauge Theories 

In this section we extend the matrix integral techniques to analyze J\f = 1 gauge 
theory with SO(N) and Sp(N) gauge groups and adjoint matter; which we first pub- 
lished in llACH + 03l . By a careful consideration of the planar and leading non-planar 
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corrections to the large M SO(M) and Sp(M) matrix models, we derive the matrix 
model free energy. We do this both by applying the technology of higher-genus 
loop equations of [ACKM93, Ake96| and by straightforward diagrammatics (see e.g. 
HBIPZ78ll£icS5K 

As for SU(M), we find that the loop equation for the resolvent of the matrix model 
(which is a Dyson-Schwinger equation for the matrix model correlation functions) 
describes a Riemann surface which is identified with a factorization of the spectral curve 
of the J\f = 2 gauge theory. 

In section 15.3.21 we discuss the application of the higher-genus loop equations to 
the computation of the MP 2 contribution to the free energy. The loop equations take the 
form of integral equations which give recursion relations between the contributions to 
the resolvent at each genus. They suggest a very simple result for the MP 2 contribution 
in terms of the sphere contribution. We verify this relationship by explicitly enumerating 
ribbon diagrams with several types of vertex; as in the previous section, the combina- 
torics of these additional diagrams combine to reproduce the expected physical result, 
we find that the contribution to the free energy T\ from IRP 2 and JF from 5* 2 are related 
by 



Tx = ±q^, (5.37) 
where S is half of the 't Hooft coupling for the SO / Sp component of the matrix group. 



We determine the proportionality constant q from the diagrammatics to be q 



4 - 



Our results suggest a refinement of the proposal of Dijkgraaf and Vafa [DV02c| for 
the effective superpotential in the case of SO and Sp gauge groups. We find that 

dJ- 

W eff = Q m -^- + Q Q5 g + TS } (5.38) 
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where Qds is the total charge of D5-branes, Q05 is the total charge of 05-planes, JF 
is the contribution to the matrix model free energy from diagrams with the topology 
of a sphere and Q is proportional to T\, the contribution to the free energy from RP 2 
diagrams. We use (15.381) to obtain results consistent with gauge theory expectations. In 
particular, the subleading correction to the matrix model restores consistency with the 
requirement that there is a degeneracy of the massive vacua of the gauge theory given 
by h, the dual Coxeter number of the gauge group. 

5.3.1 The classical loop equation 

We first consider the saddle point evaluation of the one matrix integral for SO(M) or 
Sp(M) matrices. Our discussion is analogous to that of section 15.11 and consists of 
obtaining a loop equation for the resolvent. In the next section, we will formulate a 
systematic method for obtaining the g s corrections to the classical solution. 



The partition function for the model with one matrix $ in the adjoint representation 
of the Lie algebra of G = SO(M) or Sp(M) is 



In Appendix lA.il we collect results that are useful for SO / Sp groups, but here we shall 
discuss only the SO (2M) group in detail. 

In the eigenvalue basis, the integral over an SO(2M) matrix is given by 




(5.39) 




(5.40) 
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In terms of the number of eigenvalues Mj in the neighbourhood of the critical point Xi, 
define the 't Hooft couplings 



So = 9s^y, S i = g a M i . (5.41) 
The effective action for the gas of eigenvalues is given by 

i<j gs i 

Note that W is now a polynomial of order 2n with only even powers; this is because the 
trace of an antisymmetric matrix vanishes 5 . 

This action gives rise to the classical equations of motion 

E tA^ - 7 W ' {K) = °- (5 - 43) 

Defining the resolvent 

i 1 

allows us to rewrite the equations of motion as: 

uo{x) 2 - g s - u/ (x)j + f{x) - 2cu (x)W'(x) = 0, (5.45) 



where 

is a polynomial of order 2n — 2 with only even powers, i.e., it has n coefficients 



ft \ s^ 1\W\\)-2xW\x) 



5 In principle W($) could also contain the Pfaffian, but we will omit this case. 
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In the small g s limit, (15.451) reduces to 

u {x) 2 + f(x) - 2u {x)W'{x) = 0, (5.47) 
which may again be written in the form 

y 2 - W'(x) 2 + fix) = 0, (5.48) 

via the change of variables 

y{x) = lu (x) - W\x). (5.49) 

The force equation is now 

2y{\) = -g s ^ (5.50) 

where the factor of 2 comes from the fact that the force is acting on an eigenvalue and 
its image. 

In terms of the eigenvalue density p(A), 



l x 2 — A 2 J \x — X x + X J J_ 00 x — A 

the filling fractions are given by 

So = -j— I y(x)dx, 

4« J An 

x r (5-52) 



Si = — / y(x)dx , % > 



Note that we only integrate around half of the cycle A because of the orientifold pro- 
jection. 
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a) 



b) 



Figure 5.2: Feynman rules for the SO and Sp matrix models: a) untwisted and b) 
twisted propagators 

For SO(2M+ 1) and Sp(M), one can easily see that JF and the Riemann surface are 
the same as in the case of SO(2M). These gauge groups are distinguished at subleading 
order in the g s expansion of the free energy; in the next section we will determine the 
leading contribution to the free energy from unoriented diagrams. 

5.3.2 g s corrections and loop equations 

We can now expand the free energy in terms of ribbon graphs as before. The propagator 
of the SO(M) matrix model is 



Thus, the ribbon graphs now have the possibility of "twisted propagators" as well as the 
previous untwisted propagators (see figure l5T2b : an important point is that the twisted 
propagators comes with a relative minus sign. The twisted propagators can give rise 
to non-orientable ribbon graphs, so the topological expansion includes a sum over dia- 
grams that may be embedded in non-orientable Riemann surfaces. As before, the overall 
power of g s associated to a ribbon diagram is 



where g denotes the genus, and c denotes the number of cross-caps of the Riemann 
surface on which the diagram is enscribed. 



(5.53) 



9. 



-X - „2g+e-2 



(5.54) 
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5.3.3 The resolvent 

We shall now review the general technique of loop equations [ACKM93 , Ake96|, which 
is an iterative procedure to calculate corrections to the partition function of the higher 
order in g s . Central to this procedure is the loop operator, defined as 

^ (x)= -gJMi-- (5 - 55) 

The resolvent, which is the generating functional for the single trace correlation func- 
tions of the matrix model, is defined as 

^H^Tr^Wf;^ (5-56) 

k=0 



Using the identity 



the resolvent can expressed as 



-{2k)^-F = g s (Ir<$> 2k ) ) (5.57) 
ugh 



u (x) = -£-(x)F+-, (5.58) 
av x 



where we used S = Si = g s M. We are using the variables g s and S since we are 
working in the small g s limit with S fixed. As mentioned before, the perturbative part 
of the free energy has an expansion in g s of the form 

fpert. = J2g 2 s s+c ~ 2 F 9 , c (5.59) 

We will be interested in calculating the first two terms in this expansion, which are the 
contributions from diagrams with the topology of S 2 and MP 2 , although the analysis 
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can in principle be extended to all orders to study gravitational corrections to the gauge 
theory superpotential. The resolvent has a similar expansion 



W^ Ss 2 %(4 (5.60) 



The asymptotic behavior at infinity of the oo 9yC is clear from the definition of uo(x) 

S 



u ,o( x ) =- + 0(x- 2 ), 

X 

uj gyC (x) =0(x- 2 ), 2g + c>0. 

Using this fact and the existence of the genus expansion, we can write 

d S 

Wg,c( X ) = ^/^ :F 9^ 2 9 + C>0. 



(5.61) 



(5.62) 



These equations determine the dependence of T gfi on the coupling constants. There is 
still an additive constant that is undetermined, but this is physically meaningless. In 
the next section we will derive the loop equation, which will provide us with recursion 
relations to calculate uo 9tC as functions of the coupling constants gj appearing in the 
matrix potential W(M) = J2j yTrM- 7 . For the rest of the discussion, we denote uj ,o 
by uoq and a; ,i by uo\. 

5.3.4 The loop equation 

In this section we will derive an important recursion relation between the different per- 
turbative contributions uo 9yC to the resolvent. The loop equation can be derived by per- 
forming an infinitesimal reparametrization of the matrices $ in the matrix integral and 
using the fact that the integral is trivially invariant under reparametrization of $. Let us 
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reparametrize $ by 



$ = = $'_ e y£J__ ( 5. 6 3) 

\ X * / odd fc=0 X 

<i* = *- £ ££^_ (5.64) 

fc=0 1=0 

where we only take the odd/even powers of $' in order to preserve the SO/Sp Lie 
algebra. The Jacobian for this reparametrization, keeping only lowest powers of e, is 
then 

e / 1 \ e 1 



J($') = 1 - - Tr — + — Tr — . (5.65) 



The action transforms as 



= TrJV ( $' - ( — — ) ) = TrW(&) - eTr ^ ^ \ (5.66) 



^- $ 7odd/ x-$' 

Inserting this into the matrix integral, the invariance under the small variation of $ yields 
the identity 



1 \ 2 1^ 1 

Tr Tr- 



x — 7 X x — $' 



- — TrW(*') 

e as 



(5.67) 



We can now make use of the identity 



d , , , , I (_ 1 \ 2 \ /„ 1 ' 2 



_ (lMl)= Tr— -<Tr— ) (5.68) 



(which is a rewriting of the steps leading to !5 . 1 6b to get the loop equation 



9s < Tr^f ) = \^xf - |^(x) + f ^(*M*). (5.69) 
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We can rewrite the loop equation using 




x - $ 




2ni x — x' 



dx' W'{x') 




(5.70) 



where C is a contour that encloses all the eigenvalues of $ but not x. In the large M 
limit of the matrix model, we get a continuous eigenvalue distribution for $ and all the 
eigenvalues are distributed over cuts on the real axis of the x -plane. The loop equation 
now reads 



We can now insert the g s expansions (15.601) for the resolvent and iteratively solve for the 
uj g , c . The zeroth and first order equations are 



The resolvent that solves the loop equations must satisfy (I5.61I) . which imposes con- 
straints on the end-points of the cuts in the a;-plane. 

Equation (15.731) is a linear inhomogenous integral equation for u)\. The homoge- 
neous equation is solved by a derivative of tu with respect to any parameter which 
specifies the vacuum, i.e., is independent of the coupling constants gj. In our case there 
are only the parameters Si, which specify the classical vacuum around which the matrix 
integral is expanded. 




(5.71) 




(5.72) 



(5.73) 
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5.3.5 Solution to the loop equations 

We now solve the loop equations (15.721) first for lu and then for uo x in the case of a 
polynomial potential 

n 

W(§) = — $2i - ( 5 - 74 ) 

Planar contributions 

In equation (15.721) . we deform the integration contour C to encircle infinity, and rewrite 
it as 

\u {xf = w\x)^{x) + l x / ^VH(*0 (575) 

2 Jc x x - x< 

Assuming that uj (x) has k cuts in the complex x-plane, we make the ansatz 



oj ( x ) = W'(x) - M(x) 



\ 



2k 



Y[(x- Xi ), (5.76) 



i=l 



where M (x) is an undetermined analytic function at the moment. Here the end points of 
the cuts, denoted by the Xj, are unknown and have to be determined. It is clear that if we 
have the maximum allowed number of cuts, k = 2n— 1, the function M(x) is a constant. 
The loop equation determines M in this case to be the coupling constant g n . For the 
SO/Sp models the eigenvalues come in pairs, and the total number of "independent" 
cuts is n. There is one cut [— x , x ] centered around zero, and the other cuts come in 
pairs [x 2 i-i,x 2i ] and [— x 2 j, —X2i-i\. Note that the cuts are simply the projections of 
the S 3 cycles of the Calabi-Yau geometry that engineers this gauge theory, which we 
discussed in section UA\ 
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We now demand that the resolvent ujq(x) falls off at infinity as S/x (and hence 
vanishes in the classical limit S — > 0), and thus obtain n constraints 



The most general solution to these n constraints ( 15.771) is given by 

2(n-l) 

g 2 n J] (x 2 - xf) = W'(x) 2 - f(x), (5.78) 

i=0 

where f(x) is the most general even polynomial of order 2n — 2, 

n-l 

f(x) = Y,bix 21 . (5.79) 

1=0 

Note that we have now recovered the solution to the classical loop equation that we 
obtained in section 15.3.11 We now repeat the procedure outlined there and define the 
Riemann surface £ given by 



The filling fractions Si then become period integrals of the meromorphic 1-form ydx 
over the 1 -cycle Aj of S that encircles the i th branch cut 



We can then argue that the change in the free energy due to an eigenvalue tunneling to 
infinity from the i th cut is 




(5.77) 



y 2 = W\xf-f{x). 



(5.80) 




(5.81) 




(5.82) 
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This again requires the introduction of a cutoff, which can be understood in terms of the 
combinatorics of diagrams with 1 boundary and the topology of a disc or Moebius strip. 



MP 2 contributions 

Once we have the form of the solution for u (x), we can substitute it in the loop equa- 
tion, which is now a linear inhomogenous integral equation for Ui(x), 

x' — 1 ; ; = wo{x)u) X {x) - —ou (x). (5.83) 

(j x — X L*X 

We can get a natural ansatz for uj\ from the string theory expectation that T\ should 
be a derivative with respect to So of JF , 

dJ- 

?i = Qj^, (5-84) 
where q is some constant which has to be determined. Inserting this into (15.62b . we get 



/ n d i x T ^ 2j d dF 



° r u (x)--\ (5-85) 



dSo \ x 
dujQ 2q 
dSo x 



It is easy to see that solves the homogeneous part of the loop equation. The inho- 
mogenous part of the loop equation is solved by — ^ if q = — 4. 

More generally, in the case of multi-cut solutions, we could have added any solution 
to the homogeneous loop equations. This amounts to taking 

*i = X>^> (5 - 86) 



118 



such that ^2 qi = — |. However, corrections of the form ^jP- for % > should not be 
generated since these cuts represent U(Ni) gauge physics for which there should be no 
MP 2 contribution. We will give a short perturbative discussion of this in the next section. 

5.3.6 Counting Feynman diagrams with S 2 and RIP 2 topology 

For a perturbative check of the relation 

T x = ±g^ (5.87) 

we can enumerate "ribbon" graphs in the genus expansion of the matrix model. Recall 
that the genus expansion is ordered by diagram topology, with diagrams of genus g and 
c cross-caps contributing at order g~ x = g^ 2+2g+c - The coefficient q is related to the 
relative contribution of the planar (genus 0) diagrams which dominate at large M and 
the leading -g correction coming from diagrams with topology RP 2 . 

It is known that SO(2M ) and Sp(M) matrix models are related by analytic contin- 
uation M i — ► — M (for the analogous gauge theory results see [Mkr81 , CK82l lCic82l '). 
Therefore, at even orders in the genus expansion, the contribution to the matrix model 
free energy is the same for both theories, while at odd orders the Sp(M) diagrams con- 
tribute to the free energy with an additional minus sign relative to SO(2M). This fact 
determines the sign in (I5.87I ). Recall that 

x = v-p + l (5.88) 

where v is the number of vertices in the ribbon graph, p is the number of propagators and 
/ the number of boundary loops. The Feynman rules are summarized in appendix IA.3 1 
Let us evaluate the first-order quartic diagrams in fig. 15.31 The planar diagram has the 
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Planar diagrams: 




RP diagrams: 






1 

Figure 5.3: S 2 and IRP 2 diagrams with one quartic vertex, written in terms of twisted 
and untwisted propagators and as diagrams on RP 2 to show their planarity. Propa- 
gators that pass through the cross-cap become twisted. 



value 



2x ig^f^y M3 

l\4g s \2mJ 



(5.89) 



whereas the MP 2 diagram with one twisted propagator contributes 



-4xi^f^) 2 M 2 
l!4fif s \2mJ 



(5.90) 



120 



and the MP 2 diagram with both propagators twisted contributes 

l\4g s \2mJ 
Using the fact, that S = ^M, this shows that 

at the first order. We have enumerated the Feynman diagrams to several higher orders 
and higher vertices and confirmed this relationship in those cases 6 (see Appendix 1X1 for 
some examples). 

In order to describe a multi-cut matrix model (corresponding to vacua with classi- 
cally broken gauge group), we would need to use ghosts HDGKV03I to expand around 
the classical vacuum. In this prescription, one can think of the matrix model as several 
matrix models, which are coupled by bifundamental ghosts. Only one of those matrix 
models is actually an SO(M ) / Sp(M /2) matrix model, the other matrix models are 
U(Mi) matrix models. The ghosts do not have twisted propagators, so the leading con- 
tribution from the SO(M )/Sp(M /2) matrix model is again the same as for a single 
cut model. The loop equations still hold for the multi-cut model and the calculation can 
be extended to all orders. 

5.3.7 Computation of effective superpotentials 

In this section we combine the results of the previous sections to compute the effective 
superpotential of the dual gauge theories. We will find that it is necessary to refine 
the formula originally conjectured by [DV02cJ for the unoriented string contribution 

6 This combinatorial result was previously unknown to mathematicians. 
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(5.91) 



(5.92) 



to the effective superpotential. Recall that in a vacuum with coincident eigenvalues, 
there is a non-perturbative contribution to the matrix model free energy coming from 
the logarithm of the volume of the residual gauge transformations that preserve this 
vacuum. In appendix B, following [OV02|, we have included the large M expansion of 
the logarithm of the volume of the SO/Sp groups. We find that, for SO(M) when M 
is even, the non-perturbative contribution to the free energy is 

T n "P _ 1 <r-np | 1 spnp , 



9l 



1 

+ — 

9s 



S 2 log^-S 2 fUlogvr 
m V 2 



(5.93) 



S 2nS £ JN 
~2 lo § — + T^ 1 + log ^ _ log4 ) 



+ 



with a similar expression for M odd or G = Sp(M). We see that 



• F " P = T 4^ ± 2 1 ° g2 ' (5 - 94) 



where the first — /+ sign is for £0 / respectively. This is almost the same relationship 
as we found for the perturbative contributions ( 15.921 ). but it is spoiled by the log 2 term. 
This amounts to a factor of 2 discrepancy in the volume of the gauge group 7 . 

It is the non-perturbative sector, specifically the coefficient of the S^logS 1 term, 
that determines the number of gauge theory vacua, which is a main consistency test 
of the translation between matrix model quantities and the effective superpotential of 
the gauge theory. The number of vacua of a supersymmetric gauge theory is equal 
to the dual Coxeter number h of the gauge group [Wit82, Wit98|. Therefore the total 
superpotential should lead to the conclusion that S h is single-valued. 

7 This mismatch may be related to the choice of whether or not to work in the covering group. 
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Open string physics tells us that the sphere contribution to the effective superpoten- 
tial should be proportional to Qd5, the total charge of D5-branes, while the MP 2 con- 
tribution should be proportional to Q05, the total charge of 05-planes. We can express 
this by refining the suggestion of ID V02cl : 

dJ- 

W eff = Q m + Q 05 Go + rS, (5.95) 
We assume that Q is proportional to the total RP 2 free energy, 

g = a(F? p + F?). (5.96) 
Proceeding with this result, we find that 

W eff = (f ± i) S1 °g S +r 5 + ---' ^ 

where the +/— is for SO/Sp respectively. Consistency with both the closed string 
result ( 13.91) and the gauge theory 8 requires that we must have a = =p4. This was con- 
firmed by PNO03I who gave a perturbative argument along the lines of llDGL + 03l ; it 
was found to be related to the measure on the moduli space of Schwinger parameters, a 
quantity that is intrinsic to the gauge theory. 

Note that the first subleading (non-planar) contributions combine with the leading 
(planar) contributions to give a shift in the overall coefficient. This combination with 
the leading-order contributions is quite similar to the role of planar diagrams with one 

8 Note that, after including a = =p4, the effective superpotential naively suggests that for gauge group 
Sp(N/2), S N+2 is single-valued, whereas h — N/2 + 1. The resolution to this puzzle was explained 
in |Gom02|. Namely the Dl-string wrapped on P 1 has instanton number two in Sp(N/2). Properly 
accounting for this reproduces the Z-2h chiral symmetry of the dual gauge theory. 
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boundary discussed in the previous section, which conspire to soften the UV divergence 
of the planar free energy. 
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Chapter 6 



Conclusions 

In this thesis we have studied effective superpotentials for confining gauge theories with 
Af = 1 supersymmetry, focusing on theories where an underlying Af = 2 supersym- 
metry is softly broken by a tree-level superpotential. String theory provides insight into 
the structure of the vacua of these quantum theories, and a set of geometrical tools for 
computing the effective superpotential exactly, even in strongly-coupled regimes. 

The techniques we have discussed revolve around the computation of period inte- 
grals of a meromorphic 1-form on a particular hyperelliptic curve. We studied this curve 
and found that the regularization of the divergence at infinity of the contour integrals 
requires the introduction of additional fundamental matter superfields into the gauge 
theory, which cut off the domain of integration and render the calculation finite. This is 
physically pleasing, since the Af = 2 gauge theory with Nf = 2N C massive fundamental 
hypermultiplets has vanishing /3-function at high energies, indicating that the theory has 
a nontrivial conformal fixed point and is free from short-distance singularities. 

We evaluated the period integrals explicitly for the maximally confining vacua (com- 
pletely degenerate curve), and derived an explicit expression for the superpotential of an 
U (N c ) gauge theory with < Nf < 2N C fundamental superfields of arbitrary non-zero 
mass, and arbitrary tree-level interactions of the adjoint superfield $. Extremizing this 
superpotential gives the exact vacuum superpotential and the gaugino condensate, and 
agrees with previous special cases discussed in the literature. 

The Af = 1 curve may be obtained by factorizing the Seiberg-Witten curve of the 
underlying Af = 2 theory. For the Af = 2 theory with fundamental hypermultiplets, 
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we solved the factorization problem for the case when the curve factorizes completely, 
and used this solution to verify the combinatorial form of the effective superpotential in 
Af = 1 theories with fundamental matter. 

The Af = 2 gauge theories are known to have an underlying integrable structure, and 
this partly survives the soft supersymmetry breaking to Af = 1. The existence of this 
integrable system is equivalent to the statement that the vacuum of the gauge theory is 
completely characterized by the vev of the adjoint (matrix- valued) chiral superfield $; 
this matrix is identified with the Lax matrix of the integrable system, which completely 
characterizes the integrable dynamics. We considered the theory with Nf = N c fun- 
damental hypermultiplets, and found the value of the Lax matrix of the corresponding 
integrable system in the maximally confining vacuum. This form of the Lax matrix was 
not previously known. 

The geometrical techniques involving the spectral curve, while computationally 
powerful in obtaining exact results about the confining phase of supersymmetric gauge 
theories, do not have a clear origin within the gauge theory. The bridge between the 
geometrical techniques and the physics of the supersymmetric gauge theory is provided 
by the matrix models. The reason for this correspondence is that after summing the 
gauge theory diagrams order by order, the 4-dimensional loop momenta cancel and only 
the planar combinatorics survive. 

We discussed how the Af = 1 curve - and therefore the gauge theory effective super- 
potential - emerges from the study of a particular class of matrix integral, considered as 
a zero-dimensional path integral for the eigenvalues of the matrix. Specifically, the free 
energy of the matrix model receives perturbative and non-perturbative contributions. 
The former come from the planar (and for SO and Sp gauge theories, or theories with 
fundamental matter, the leading non-planar) diagrams of the matrix integral, which are 
in 1-1 correspondence with Feynman diagrams of the gauge theory. 
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The non-perturbative contributions to the matrix integral come from the residual 
gauge transformations that exist when two or more eigenvalues populate the same crit- 
ical point of the potential. These correspond in the four-dimensional gauge theory to 
classical vacua where a subgroup of the gauge group remains unbroken in the clas- 
sical theory. Expanding the volume of the gauge groups reproduces the Veneziano- 
Yankielowicz superpotential. In the four-dimensional gauge theory this superpotential 
is generated by the strong-coupling dynamics of the gauge field, and the only existing 
derivations come from anomalous symmetry constraints. However, since the gauge field 
does not appear in the matrix integral (in a sense, it is integrated out), there is no com- 
plication from strong coupling and the contribution may be read off from the asymptotic 
expansion of the volume of the unbroken gauge group. 

The various techniques we have used to study the effective superpotentials may 
be characterized as geometrical, algebraic and combinatorial in nature. Each of them 
involves the spectral curve, but highlights a different aspect of its structure. This struc- 
ture is in turn reflected in the structure of the vacua of the J\f = 1 gauge theory. 

These techniques teach us about confinement and other non-perturbative phenomena 
in the Af = 1 gauge theories; for example in theories with an adjoint chiral superfield 
(which contains a scalar field), confinement of the low-energy gauge theory is asso- 
ciated to condensation of the magnetic monopoles of the gauge theory, and moreover 
the exact value of the monopole condensates can be calculated. Extremizing the gauge 
theory effective superpotentials gives exact non-perturbative results about the vacua of 
the theory, such as the values of the gaugino condensates associated to chiral symmetry 
breaking. We are therefore able to obtain exact results about confining theories that are 
believed to have many similar properties to non-supersymmetric Yang-Mills theory and 
QCD, for which analytical results are lacking. 
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Though string theory is fully a theory of gravity and other fundamental forces, it is 
commonly the case that the effects of gravity can be consistently decoupled, and string 
theoretical techniques can be used to study the remaining low-energy supersymmetric 
particle interactions in isolation. Thus, if supersymmetry is realized in nature at an 
experimentally accessible energy scale, then - whether or not string theory is the correct 
unified theory of quantum gravity and fundamental forces - string theory has provided 
tools that will be useful for understanding aspects of physics in supersymmetric regimes. 
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Appendix A 



Matrix Integral Measures and 



Determinants 



In this appendix we collect some results on the group measure and adjoint action which 
are needed to do computations in the matrix models. 

A.l The group measure for general matrices 

We wish to compute the Jacobian for the transformation from certain matrices $ to their 
eigenvalues. This can be derived by a group-theoretic argument. In terms of the Cartan 
generators if; and ladder operators E a , for the algebra of the group G, satisfying 



[H\ E a ] 



(A.l) 



we can diagonalize a matrix $ 




(A.2) 



We will define parameters t a so that 




(A.3) 
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The infinitesimal variation of $ can then be written as 



d$ =U ] 



(A.4) 



U. 



We now calculate the metric on the Lie algebra 



Tr d$ d& = + Yl dta dt P ( Yl Xi0ii ) ( Yl ) Tr EaE P- (A - 5) 



Using the identity 



TirETE = C{r)5 a 



+13,0 



(A.6) 



where C(r) is a representation dependent constant, we can simplify the second term in 
equation dA.51) to 



Q: 

Up to numerical factors, the Jacobian is 



\dt a \ 



(A.7) 



a(a) = n 



a>0 



(A.8) 



We list the expressions for the roots and the corresponding determinants for the 
different classical groups in Table IA.l1 



A.2 Asymptotic expansion of the gauge group volumes 

We now compute the asymptotic expansion of the volume of the gauge groups, which 
normalizes the partition function of the matrix model and provides the nonperturbative 



139 



G 

Roots 


J(A) 


A N -i 

d - ej (i ^ j) 


ll^-A,) 2 

i<j 


B N 

±e; ± ej (i ^ j), ±ti 




Cn 

j- (±e< ± ej ) (i ^ j), ±V2e, 


ito-a^eia 2 


D N 

±d ± ej (i ^ j) 


n(A 2 -A 2 ) 2 



Table A.l: The roots and the formulae for the Jacobians associated to the classical 
groups. 



contribution to the free energy. The volumes are given by WOV02I : 

, , . . ^N(2ti)I n2+1 2 N - 1 

2 N+1 (2ir) N2+N -* 
TO " 50 ' 2iV+1 >»- (2 J V-l)!(2 J V-3)!^.3!ir 



vol(SO(2A0) 
\ol{Sp(2N)) - 



^/2{2it) 



A 



(A.9) 



(2AT-3)!(2AT-5)!...3!1!(A^-1)! 

2-N( 27T ) n2 + n 



(2iV-l)!(2JV-3)!...3!l!" 



We are interested in the large N asymptotic expansion of the logarithm of the vol- 
umes in order to compute the non-perturbative contribution to the free energy. Following 
IIOV02I . we introduce the Barnes function 

G 2 {z + 1) = T{z)G 2 {z), G 2 (l) = 1. (A.10) 
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Using the doubling formula for T(z), 

T(2z) = 2 2 ^ 1 vr-^r(z)r(^ + i), (A.ll) 
and (IA.10I) . can evaluate the denominator of the volume factors 

G d (N) = (2N - 1)! . . . 3!1! = —±—2 N ( N+1 '>G 2 (N + l)G 2 (iV + 1) (A.12) 

Using the Binet integral formula 

1 1 f°° tan(-) 

logr(z) = (z- -)log* - z + - log 2tt + 2 J -^rf^dt, (A.13) 

the asymptotic expansion of G 2 (n) is 

jy2 i 3 1 

log G 2 (N + 1) = — log AT - — log N - -N 2 + -iVlog27r + 0(1). (A. 14) 

I J. _ 

By expanding log(iV — a) for large N, we obtain 



log G d (N) =N 2 logiV + iV 2 (-^ + log 2) 

1 1 N 

+ -iVlogiV - — logiV + - (log4vr - 1) + O(l). 



(A. 15) 
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Putting all of this together, we find that 

log \ol(SU(N)) 

= —N 2 logiV + — logN 
12 

+ ^ 2 + iiV 2 log27r + 0(l), 
log vol(50(2AT + 1)) 



= —N 2 logN + N 2 {- + log7r) 
2 



1 IN 
- -N log N + — log N + -(1 + log 4 + log tt) + 0(1), 

^ ^4 ^ (A16) 



logvol(50(2iV)) 



= -iV 2 log AT + N 2 {- + log7r) 
2 

1 1 N 

+ -AHogiV + — logiV + -(-1 + log4 - logvr) + 0(1), 

log vol(Sp(2iV)) 

= -iV 2 log + iV 2 (- + logvr) 
2 

1 1 TV 

- -NlogN + — log TV + -(1 - log 4 + logTr) + O(l). 

A.3 Matrix model Feynman rules and enumeration of 
diagrams 

We want to perturbatively evaluate the matrix integral 

J d^e^ Trw ^\ (A.17) 

where the potential W is given by 

oo 

W($) = ^2^$ 2j (A.18) 

3=1 23 
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and $ is a real antisymmetric M x M matrix. We can write this as 



/ rf$ exp -Tr ( -$ 2 + 



(A. 19) 



where m — g\. Expanding the exponential leads to traces of integrals of the form 



1 Trmj.2 

dSe* 11 ** $ mini • • • $ mfe „ fe 



d 



d 



(9 7 (9 7 

x ~>' J m\ri\ <^' J m k n k 



d$ exp 



9s 2 2 



j=o 



(A.20) 



This integral can now be evaluated, leading to 



iU(M-l) 

27rg s \ 2 d 



d 



m 



(9 7 (9 7 

L/ <-'mi?ii x ~>' J m, k n k 



__2s.Tr/ 2 



) J=0 



(A.21) 



Differentiating step by step gives rise to expressions like 



'mn 



7 t p-|^Trj 2N \ 



9J m „ V2m" mini 2m ™ 



_/ r r _r r \ _^£_ T J__7 -|^Trj 2 

2m 2m 2m 

+ ■■■ 

9s T 9s T 9s ,, r _x x \ -f^Trj 2 

2m ' ' ' 2m 2m °mn k Onm k )C m 

,_9^t_9^j _9^ T P -#^Trj 2 

2m 2m 2m 



(A.22) 



The indices m^ and n» are contracted in traces as given in the interaction which can be 
interpreted as forming vertices. The combinatorics can then be interpreted diagrammat- 
ically; one must connect all the legs of the vertices in all possible ways with untwisted 
and twisted propagators. Each twisted propagator contributes a factor of (—1). 
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The rules for evaluating a diagram are then: 



Each kind of vertex with multiplicity Vj contributes a factor of yi{-^-) Vl • 



Each propagator contributes a factor of 



Each twisted propagator contributes an additional factor of (— 1). 



Each index loop contributes a factor of M = —. 



The combinatorial factor of a diagram can be computed by counting all topologically 
equivalent ways in which the legs of the vertices can be connected. This has some sub- 
tleties, since some diagrams with twisted propagators can actually be planar. To handle 
this, we make use of the technique described in HCic82H to draw unoriented diagrams 
(see also HMW031 IMY02I for recent work on non-orientable ribbon diagrams in the 
mathematical literature). 

An WF 2 can be drawn in the plane as a disc, where antipodal points on the boundary 
are identified. WF 2 diagrams can then be drawn on that disc with some propagators 
going through the cross-cap at the boundary. The propagators going through the cross- 
cap are twisted propagators, whereas all the others are untwisted propagators. 

We can now also draw a planar diagram on the MP 2 . If it has more than one vertex, 
we can push one or several vertices through the cross-cap without destroying the pla- 
narity, but all the propagators going through the cross-cap are now twisted propagators. 
This operation contributes a multiplicative factor of 2 1 " 1 to the number of planar dia- 
grams at each order v. See Figure 1531 for the enumeration of diagrams with 1 quartic 
vertex. 

Using the relation between p and the number of vertices Vi of valency i 




(A.23) 
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the contribution of planar diagrams to the free energy of the SU(M) matrix model is 
given by 




(A.24) 



where the sum is over diagrams with v vertices of valence 2n, d v is the number of 
planar diagrams at each order, and / counts the number of boundary loops of the ribbon 
graph. The propagator for SU(M) theories is twice that of the SOISp theories. In the 
second line we have simplified using (15.881) and (IA.23I) . The number of diagrams of 
topology S 2 (i.e. planar diagrams) in SU(M) matrix theory with a quartic potential is 
given by IBI PZ78ll 

4 4) = ^7— = 2, 36, 1728, 145152, .... (A.25) 
(y + 2)! 

We are not aware of explicit generating functions for other vertex valences 2n, but these 
diagrams can be enumerated by computer to the desired order. 

If we now include twisted propagators (i.e. enumerate planar diagrams in the SO 
or Sp matrix models), there is an extra contribution to the set of planar diagrams com- 
ing from vertices that have been "flipped", converting untwisted to twisted propagators 
according to the rule described above. 

00 ,(n) 00 ,(n) 

F Q = y ^L_(^W^)p M ' = Y (a.26) 
/ — i v\ nq s 2m v\ nq s 2m 

v=l v=l 

, m 1 (2v - l)\24: v 

" = 2 (^ + 2)! = 2 ' ' 6912 ' 1161216 ' ■ • • • (A " 27) 
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A similar expression exists for the IRP 2 free energy 



oo An) 



7-1 



oo An) 

y ^_ ( ^ r( ^ ) > M i-(i-f )v (A 2g) 



Here the number of diagrams d^' is counted with a minus sign for each twisted prop- 
agator 1 . The relevant planar and RP 2 diagrams were enumerated by computer up to 
4 vertices with a quartic potential Wf ree ~ $ 4 , to 2 vertices with a sextic potential 
Wt ree ~ $ 6 , and for a single vertex with a potential of degree up to 16. The results are 
summarized in Table IA.21 and verify the desired relation: 



'Gaussian Ensembles are matrix models that have been well-studied in the physics and mathematics 
literature. The Gaussian Orthogonal and Gaussian Symplectic Ensembles also contain non-oriented rib- 
bon diagrams with twisted propagators, however the propagator is (T b a T£) ~ S ac Sbd + 5ad5bc, '•<?., there is 
no relative minus sign between the two terms. This corresponds to counting RP 2 diagrams with a positive 
sign always. Therefore the free energy of the Gaussian Ensembles differs from that of the Lie Algebra 
matrix models at sub-leading orders in the genus expansion. 



ldT< 



(A.29) 



2dM' 
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Diagrams with quartic vertices: 



Gauge group 


Topology 


v = 1 


v = 2 


v = 3 


w = 4 


SU 


S 2 


2M 3 


36M 4 


1728M 5 


145152M 6 


SO/Sp 


s 2 


2M 3 


72M 4 


6912M 5 


1161216M 6 


SOISp 


MP 2 


-3M 2 


-144M 3 


-17280M 4 


-3483648M 5 



Diagrams with sextic vertices: 



Gauge group 


Topology 


v = 1 


v = 2 


SU 


S 2 


5M 4 


600M 5 


SOISp 


s 2 


5M 4 


1200M 6 


SOISp 


MP 2 


-10M 3 


-3600M 5 



Table A. 2: Contribution to the free energy of the SU/SO/Sp matrix models at planar and 
MP 2 level, for quartic and sextic potentials. The first few terms in the perturbative expansion 
are listed, corresponding to the number of diagrams with increasing number of vertices 
(equivalently loops). 
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Appendix B 

Emergency Proof Techniques for 
Physicists 

1 . Proof by Intimidation 

Best applied to Graduate Students. 

2. Proof by Divine Revelation 
See llWitl . 

3. Proof by Exhaustion 

Keep going until your entire audience has fallen asleep or lost interest. 

4. Proof by Vigorous Gesticulation 

If you think about it, this one is pretty similar to (flTl) . 

5. Proof by Extrapolation 

Prove the result in a certain limit, then assume it holds true over the entire param- 
eter space. 

6. Proof by Physicality 

Any result you dislike is declared to be unphysical and cast out. 
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7. Proof by Approximation 

Keep approximating until the result becomes trivial, then go back and fill in some 
of the gaps. 

8. Proof by Assertion 
"It is clear that..." 

9. Proof by Dimensional Transmutation 

c = h — a! — 1 

10. Proof by Bastardized Notation 

It's easier than getting the mathematics correct. 

11. Proof by Obscure Citation 
For full details, see [Mor56|. 

12. Proof by Omission 
"It can be shown that..." 

13. Proof by Peer Pressure 

"It should be completely obvious to every reader that.." 

14. Proof by Recursive Citation 

Instead of citing the proof, cite a paper which refers to the proof, and iterate. 
Bonus points if you can introduce a cycle into the graph of references. 
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15. Proof by Conclusion 

The consequences of the result are so profound that it must be true. 

16. Proof by Redefinition 

Derive a result which is manifestly true, then redefine the meaning of the symbols 
and continue to use the result. 

17. Proof by Analogous Reasoning 

Compare the situation to a different, but vastly simpler one for which the result is 
true, and argue that the general case should have similar properties. 

18. Proof by Trivial Limit 

The result reduces to the correct one in a suitably nice limit. 

19. Proof by Rational Approximation 

2 = 7T = I = -1 = 1 

20. Proof by Opressive Citation 

Cite an unrelated 100-page paper on the assumption that no-one will search 
through the entire thing for the proof. 

21. Proof by Complication 

Spend 98% of the paper deriving impenetrable technical results, and tie together 
at least 10 different threads which miraculously produce your result on the last 
page. 
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22. Proof by Notational Gymnastics 

Change your notational conventions at least three times to distract any hopes of 
pursuit by the reader. 

23. Proof by Example 

Show the n = 1 case. 

24. Proof by Negative Reasoning 
The opposite of the result is false. 

25. Proof by Forward Citation 

"We intend to present a proof of this result in work which is currently under prepa- 
ration." 

26. Proof by Intuitive Diagram 

Draw a pretty enough picture and you can prove anything. 

27. Proof by Deception 
Watch the hand... 

28. Proof by Numerology 

We get the same numbers from two unrelated computations, so they must have the 
same meaning. 

29. Proof by Profanity 

Should probably only be used as a last resort. 
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30. Proof by Association 

Tie the desired result to an unrelated discussion of obviously true material. 

31. Proof by Universal Convergence 

It gives the right answer, so we'll let the mathematicians figure out why. 
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